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I. INTRODUCTION 

THREED is the code name of a computer program which has been developed 
for performing absolute orientation by the method of three-dimensional projec- 
tive transformation. It has the capability of performing complete error analy- 
sis on the computed transformation parameters as well as the transformed 
coordinates . 

The accuracy of absolute orientation depends on the following factors: 

1) accuracy of the model coordinates, 

2) accuracy of the ground controls, 

3) density and distribution of control points, 

4) size of the area, and 

5) scale of the stereo model. 

Program THREED was coded in FORTRAN IV computer language for the IBM 
System 360/75 computer at the University of Illinois at Urbana-Champaign. It 
may be used to perform any one of the following functions: 

1. To perform absolute orientation. 

The program takes as input the model coordinates of a set of model 
points and the ground coordinates of a group of control points. Both 
the model and the ground coordinates of the control points can be weighted 
individually according to their variance-covariance matrices. The pro- 
gram computes the seven transformation parameters (Xy, Yy, Zy, w, 41, k 
and scale) and their estimated standard errors. The program also trans- 
forms the model coordinates of any pass point into the ground reference 
system and determines the standard errors of the transformed coordinates. 

2. To study accuracy of absolute orientation by simulation. 

In simulation application, the program takes as input 1) the 
ground coordinates of a set of control points; 2) the variances of 
the ground control coordinates; 3) the variances of the model coor- 
dinates; and 4) the scale of the model. The program then generates a 
set of model coordinates for the given ground points and perturbs 
them according to the specified accuracy of the model points. It then 
performs a regular absolute orientation solution and outputs the estimated 
standard errors of the seven absolute orientation parameters. 
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3. To determine the uncertainty in the orientation of the surface 
defined by a set of triangulated pass points. 

The direct output of a phototriangulation solution is the ground 
coordinates of a set of pass points and their standard errors. Program 
THREED can be used to determine the uncertainty in the orientation 
(X-pj Y-p Zp id, 0, k and scale) of the surface defined by the set of 
pass points. The program takes as input the pass point coordinates and 
their standard errors. It then generates fictitious ground control 
coordinates which are translated, rotated and may have different scale 
with respect to the pass point coordinates. The pass point coordinates 
are then treated as model coordinates and transformed into the ground 
system. By assigning very small standard errors to the fictitious 
ground coorindates, the standard errors of the computed transformation 
parameters then reflect the uncertainly in the orientation of the 
surface defined by the set of pass points. 

Program THREED was developed for the purpose of a research study on 
the treatment of control data in lunar phototriangulation. The project 
was sponsored by NASA - Lyndon B. Johnson Space Center under Contract No. 

NAS 9-12446. The application of program THREED in studying the accuracy the 
lunar phototriangulation was reported in the final technical report of this 
project (1). Nick G. Yacoumelos, presently an Assistant Professor at Lowell 
Institute of Technology, was the research assistant on this project and was 
responsible for the coding and testing of program THREED. 

2. MATHEMATICAL FORMULATION 

The program THREED is based on the equations for three-dimensional 
projective transformation which are as follows: 


(1) Wong, K. W. , "Treatment of Control Data in Lunar Phototriangulation," 
Civil Engineering Studies, Photogrammetry Series No. 37, University of 
Illinois at Urbana-Champaign, Urbana, Illinois, 61801, January 1974. 
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where 



X 


model coordinates of point j 
ground coordinates of point j 
three translations 
scale factor 


m-|i = cos $ cos < 

m-|2 = cos to sin k + sin to sin <t> cos k 

m^ = sin w sin k - cos to sin 4> cos < 

m^i = -cos (j> sin k 

ni22 = cos to cos k - sin to sin <#> sin ic 

TTI23 = sin to cos < + cos to sin 4> sin < 

m 31 = sin <f> 

m 32 = -sin to cos 4 1 

m 33 = cos to cos (J) 

After linearization by first-order approximation, Eq. (2.1) may be expressed 

as follows: 



i.e. Vj ♦ B j4 ♦ BjAj = 


( 2 . 2 ) 
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The model coordinates of each control point generates one set of equations 
as in Eq . (2.2). For m control points, the complete set of observation equations 
may be expressed as follows: 



i.e. V + Ba + Ba = e (2.3) 

In order to permit flexible weighting of the seven transformation 
parameters as well as the ground control coordinates, one set of observation 
equations is introduced for each. The set of observation equations for the 
transformation parameters are as follows: 
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where X*, yS, ... and x° are approximated parameters; and xS°, yS°, ... 
and x are measured parameters. In matrix notation, this set of equations 
may be simply written as follows: 

V - A = £ 

(7,1) (7,1) (7,1) 


(2.4) 
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The observation equations for the jth ground control point are as follows: 
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y 0 yOO 
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Again, the superscript (o) denotes approximation parameters and the superscript 
(oo) denotes measured parameters. Equation (4.5) may be simply written as 



The complete set of observation equations for all control points are as 
fo 1 1 ows : 
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i-e. V - A - e (2.6) 

Combining Eqs. (2.3), (2.4) and (2.6) yields the following observation 
model : 
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i.e. V + Ba = C (2.7) 

The normal equation is then as follows: 

(B T WB)a = B T WC (2.8) 

where W is the weight matrix for the observations. An iterative solution 
procedure must be followed. An initial set of approximate values is assigned 
to all unknown transformations parameters. The solution solved for the 
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corrections and then apply the corrections to the approximations. The 
solution is iterated until a stable solution is reached. 

After the last iteration, the variance-covariance matrix (a T ) of 
the computed transformation parameters is computed by the following ex- 
pression: 

a T = a 2 (B T WB) _1 (2.9) 


where a Q is the variance of unit weight. The ground coordinates of all 
other model points and the corresponding standard errors are computed by the 
following expressions: 

Y j " T * m 12 x j + m 22 y j + ™32 Z j* + Y T 
z j = I (m i3 x j + "a* + "W + z t 

°x . 2 - < m n x j + Wj + "si 2 /^. 2 

+ [-(sin<j>cosK)x. + (sin^sinicjy . + cos<f>z .] 2 x~ 2 a 2 
J J J <t> 

+ [-cos^sinicx. - cos^cosicy,] 2 ^" 2 ^ 2 

J J K 

2-2 2 2 -2 2 2 -2 2 2 

"’ll 1 "X. + *Si x °y. +B 3i x °z. + % 

J J J T 

Y 2 ' <»12 x j + Vi + m 32 z j) 2x ' 4o x. 2 
J j 

+ [{-sintosinK + cosa)Sin 4 >coshc)x . + (-sinwcosK - cosa)Sin4>siriK)y . 

J j 

2-2 2 

- COS(0COS<J>z .] \ a 


2-2 2 


+ [(sinwcos<{>cosK)x. - (slncocos^sinic)y . + sinu)Sin<f>z .] a a 

J j j 
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2-2 2 

+ [(cosujcosk: - sinwsin<f>sinK)x. - (cosajsin< + sinusitHcosiOy.;] A a 

J J K 

■ 2 -2 2,2 -2 2 , 2 -2^ 2,2 

+ m 12 x °X. + m 22 X °Y - + m 32 X Z. + Y T 
J J J 

% 2 = < m 13 x j + “2* + ra 33 z / X 'V 
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+ [(cosojstriK + sinwsin<f>cosK)x. + (cosojcosk - sinwsinijisinKjy . 

J J 
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- simoCOStjjZ.] X a 
J 0) 

2-2 2 

+ [(rCOSu)COSij)COS<)x j . + (cosucos^sinKjyj - cosusin^Zj] x 
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+ m^ 3 x“ a x + m 23 x_ a Y + m 33 x a Z + a Z (2.10) 

J J 3 

The m.. terms in the above equations are defined as in Eq. (2.1). 

• J 


3. INSTRUCTION FOR DATA INPUT 


The first three cards in the data deck are the same for all three 
applications of the program THREED. These cards are specified as follows: 

Card 1. Parameter CArd (6110) 

Col. 

1-10 JJ To specify type of problem 

= 1 to study accuracy of absolute orientation by 
simulation 

= 2 to perform absolute orientation 
= 3 to determine uncertainty in the orientation of a 

surface defined by a set of triangulated pass points. 
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11-20 NP 


21-30 NCP 

31-40 MAX IT R 
• 41-50 JPRINT 


51-60 I PRINT 


Card 2. (F20.10) 

Col . 

I- 20 SIGO 

Card 3. Variances of 
Col . 

1-10 

II - 20 

21-30 

31-40 

41-50 
51 -60 
61-70 



Total number of points in the model, including both 
control points and points to be transformed into the 
ground reference system. 

Number of control points, i.e. points for which the 
model and ground coordinates are both known. 

Maximum number of iterations allowed. 

= 0 Full variance-covariance matrix to be input for 
ground coordinates. 

= 1 Diagonal elements of variance-covariance matrix to 
be input for ground coordinates 
= 0 Full variance-covariance matrix to be input for 
model coordinates 

= 1 Diagonal elements of variance-covariance matrix 
to be input for model coordinates 


Standard error of unit weight 

Input Transformation Parameters (7E10.3) 

Variance of input scale 
Variance of X-translation 

Variance of Y-translation 

Variance of Z-translation 

Variance of w-rotation 
Variance of ^-rotation 
Variance of x-rotation 



-9- 


The contents of the remaining input cards will depend on the purpose 
for which program THREED is to be used. The input data format for the three 
different applications of the program will be described separately on the 
following sections. 

3.1 To Perform Absolute Orientation (JJ = 2 in Card 1) 

Subroutine INPUT2 governs the data input for this application. 

Card 4. ID numbers of the four corner control points (4110) 

Col. 

I- 10 N1 A control point located at the upper left-hand corner of 

the area 

II- 20 N2 A control point at the upper right-hand corner 

21-30 N3 A control point at the lower right-hand corner 

31-40 N4 A control point at the lower left-hand corner 

These control points will be used by the program to compute preliminary 
approximations to the transformation parameters. 

Card 5 Sequence. Ground Coordinates of the Control Points (I5,3F15.3) 

There should be one card for each control point giving a total of 
NCP cards in this sequence, with NCP being specified in Card 1 of the deck. 

Col. 

1-5 ID ID number of control point 

6-20 X X-coordinate 

21-35 Y Y-coordinate 

36-50 Z Z-coordinate 

Card 6 Sequence. Model Coordinates of All Points (I5,F16.3,F14.3,F15.3,28X,I2) 
There should be one card for each model point giving a total of NP 
cards in this sequence. 
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Col. 

1-5 ID ID number of model point 
6-21 x x-coordi.nate 
22-35 y y-coordinate 
36-50 z z-coordinate 

79-80 II =1 if this is the last card in the sequence. 

Card 7 Sequence. Variance-covariance Matrix of the Ground Coordinates 

(I5,9E8.2,I3) 

There should be one card for each ground control point giving a total 
of NCP cards in this sequence. 

Col. 

1-5 ID ID number of control point 

6-13 a x 2 

14-21 o XY 

22-29 a xz 

30-37 ct X y 

38-45 o Y 2 

46-53 a YZ 

54-61 ct xz 

62-69 cr YZ 

70-78 a z 2 

79-80 II =1 for last card in the sequence 

Card 8 Sequence. Variance-Covariance Matrix of Model Coordiantes (I5,9E8.2,I3) 
There should be one card for each model point, and a total of NP cards. 
Same format description as in the Card 7 Sequence. 

3.2 To Study Accuracy of Absolute Orientation by Simulation (JJ = 1 in Card 1] 
Subroutine INPUT! governs the data input for this application. 
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Card 4. Parameter Card 
(F10.5,I10) 

Col . 

I- 10 GSCAL Scale of model coordinates with respect to ground 

coordinates 

II- 20 IAREXP Integer power of 10 which multiplies input ground 

coordinates to give the desired dimension on the 

ground coordinates for the purpose of simulation. 

The factor IAREXP provides a means of varying the dimension of the 

area without having to change the input ground coordinates. For example, 

by letting IAREXP - 2, the program multiplies all input ground coordinates 

2 

in Card 5 sequence below by a factor of 10 . 

Card 5 Sequence. Ground Coordinates of Data Points (3F15.5) 

Each card will define the approximate location of one point. There 
should be NP cards in this sequence, with NP being specified in Card 1. 
Col. 

I- 15 X X-coordinate 

16-30 Y Y-coordinate 

31-45 Z Z-coordinate 

Card 6. Variances of the Model Coordinates (3E10.3) 

Col. 

'-’0 ° x 2 

II - 20 a 2 
21-30 <r z 2 

These variances will be used to compute the weights for the model 
coordinates (x, y, z) . 

Card 7. Variances of the Ground Control Coordinates (3E1Q.3) 

Col . 

I- 10 a x 2 

I I - 20 a y 2 

21-30 o z 2 
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These variances will be used to compute the weights of the ground 
control coordinates. 

Card 8. Translations and Rotations of Model Coordinates with Resepct to 
Ground Coordinates (10X,6F10.5) 


Col . 
11-20 

XC 

X-translation 

21-30 

YC 

Y-translation 

31-40 

ZC 

Z-translation 

41-50 

0ME 

w-rotation 

51-60 

PHI 

^-rotation 

61-70 

CAPA 

■e-rotation 

The ground coordinates given in Card 5 sequence above will be translated 


and rotated according to these factors to generate the fictitious model 
coordinates (x, y, z) . 

Card 9. Seeds and Standard Deviations for the Random Realignment of the 
Ground Coordinates 2(110, F15. 5) 


Col. 



1-10 

1X1 

A seed number for generating random numbers. It should 
be an odd integer with up to 9 digits. 

11-25 

SI 

Standard deviation of the random perturbation to be 
applied to the X and Y coordinates of the ground controls. 

26-35 

1X2 

A second seed number 

36-50 

S2 

Standard deviation of the random perturbation to be 
applied to the Z-coordinates of the ground controls. 


The ground coordinates in Card 5 sequence are perturbed by the above 
standard deviations to derive the ground control coordinates for the simulation. 
These perturbations serve to create a slightly irregular pattern to the distri- 
bution of the ground controls.. 
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Card 10. Seeds for the Perturbation of Model and Ground Coordinates (6110) 

Col. 

I- 10 1X3 

II- 20 1X4 

21-30 1X5 

31-40 1X6 

41-50 1X7 

51-60 1X8 

These seed numbers will be used to generate random perturbations to 
the model and ground coordinates in the simulation. These perturbations 
represent random errors in the measurements of the model and ground coordinates. 
They will follow a normal distribution with the variances defined in cards 
6 and 7 above. 

Card 11. Perturbations of the Seven Transformation Parameters (7F10.6) 


Col. 

1-10 

DSCAL 

Scale perturbation 

11-20 

DXC 

X-translation perturbation 

21-30 

DYC 

Y-translation perturbation 

31-40 

DZC 

Z-translation perturbation 

41 -50 

DOME 

w-rotation perturbation 

51-60 

DPHI 

^-rotation perturbation 

61-70 

DCAPA 

K-rotation perturbation 


These perturbations are applied to the true values of the transformation 
parameters to derive realistic initial approximations to these parameters. 

3.3 To Determine the Uncertainty in the Orientation of the Surface Defined 
By a Set of Triangulated Pass Points. (JJ = 3 in Card 1) 

Card 4. Four seed numbers (4110) 

Col . 

I- 10 1X1 

II- 20 1X2 

21-30 1X3 

31-40 1X4 

All seed numbers must be odd integers with up to 9 digits. 
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Card 5 Sequence. Pass point coordinates from triangulation solution 

(I10»3F15.2) 

There should be one card for each pass point, and a total of NP 
(=NCP) cards with NP being specified in Card 1. 


Col. 

1-10 

ID 

ID number of pass point 

11-25 

X 

x-coordinate 

26-40 

y 

y-coordinate 

41-55 

z 

z-coordinate 


Card 6 Sequence. Standard errors of pass point coordinates (I10,3F20.10,8X,I2) 

There should be one card for each of the pass points in the card 5 sequence 

above. 

Col. 

I- 10 ID ID number of pass point 

II- 30 a x 

31-50 a 
51-70 a 2 

79-80 II s 1 to indicate last card in this sequence 


Card 7. Mean and standard deviation for the rotation parameters (2F10.5) 

Col. 

I- 10 DEG1 Mean rotations and k) in radians to be applied 

to the fictitious ground coordinates 

II- 20 DEG2 Standard deviation within which a random perturbation 

is to be generated for each of the w, $ and k rotations. 
The parameters DEG1 and DEG2 define the distribution from which fictitious 
values are generated for the three rotations to, $ and k. A set of fictitious 
ground control coordinates will be generated to have such rotations with 
respect to the pass point coordinates in card sequence 5. 

Card 8. Mean and standard deviation for the translation parameters. (2F10.5) 

Col . 

I- 10 ATRA Mean translations in X, Y and Z 

II- 20 DTRA Standard deviation of the random perturbations 
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These parameters are used to generate fictitious translations. 

Card 9. Mean and deviation of the model scale. (2F10.5) 

Col. 

I- 10 ASCAL Mean scale 

II- 20 DSCAL Standard deviation of random perturbation. 

Card 10. Variances of the ground coordinates 3(E10.3,5X) 

Col. 

I- 10 a x 2 

II- 20 a y 2 
21-30 a z 2 

2 

These parameters should be assigned very small values, e.g. 
cr y 2 = a z 2 = 1 .000 E-10. 
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3.4 CONTROL CARDS USAGE 


3.4.1 CONTROL CARDS TO LOAD THRfED ON DISK 


/♦ID TlMEB(2,0O),LlNr5*4000#lrPE«*4O0O,RE<;i0N»3OOK 
// FXrr DUMMY 

//B Dr UNIT»DISK,V0L*SER«UlU5F*,0SN*uSFR,P46n t THREED, 
// DISP«(#CATLG)'SPACE«(TRKj (30*5*2)) 

// EXfC E0RT*rE G ION*25OK 
//FORT .SYSIN DD * 

***** SOURCE DECK ***** 


/* 

// Exec LKf,OFORT* G OFlLf* , ’U$ER.P46 77 t THPEEr(THREED)"*REGlON,GO«2SOK 

//GO, SYSIN DD * 

***** TEST DATA SET ***** 

/* 


3.4,2 CONTROL C^RDS TO USE T^PEEd from dtsk 


/* ID TImEp(1»0)#LINES*10U0# REGION. 23?k 
// EXfC PROGFOPt*PROG*tHReED*PEGiOn=?32K 
//STEPLIB 00 D$N«USER,P4677,ThpEED#DISP«SHR 
//GO. SYSIN DD * 

***** DATA set ***** 


/* 
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3.4,3 CONTROL CARDS TO UPDATE THREED ON PlSK 


/*I0 TlME«(l*00>,LlNE$«l500,RrGlON»250K 
// EXFC F0RT>REGI0N«250K 
//FORt*SYSIN DU * 

***** REVISED SUBROUTINES ***** 

/* 

// EXEC LKG0F0RT»GQriLE*"USER.P4677*THPEEn« 
//LKED.SYSIN OD * 

include SYSLHOD(THREED) 

ENTRY MAIN 
NAME THREEO(R) 

/* 

// Ex EC COMPRESS 

//SySpDS DD USN.USER t p4677tTNREED#DlSPsPLD 
/* 


3,4,4 CONTROL CARDS TO RUN THREED NlTH SOURCE DECK 


tlOOO# REGION *2 3?K 

// tXEc rORTLDGO#PARM,FORT«**NOSOURcE#NOMAp",RL G ION,GO«232K 
//FORT *SYSIN DD * 


/*I^JIME»(1*0) ? LINESi 


SOURCE DECK ***** 


/* 

//G0.SY5IN OD * 
***** DATA SFT ***** 


/* 



4. FLOW CHARTS 
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4.1 Flow chart of THREED-MAIN 



READ: JJ,NP,NCP 
MAX ITR, JPR I NT, IPRINT 

( READ:S!G0 

, 

READ VARIANCES OF 
APPROXIMATE TRANSF0 
RMATION PARAMETERS 

_~1 

FORMULATE the 

DSIG Matrix 


CALL OUTPUX TO 
PRINT OUT GENERAL 
INPUT DATA 



COMMENT 


CALL TRANSF 


END 
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4.2 Flow chart of INPUTl 
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/ PRINT TRANSFORMATION "PARAMETERS 
/ THEIR APPROXIMA TIONS 

/print LIST OF GROUND AND MODEL / 
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4.3 Flow chart of INPUT2 


START 


fREAD Nl t N2,N3 > N4 

/ READ NUMBER OF TT 
AND GROUND COORDINATES 
OF CONTROL POINTS 

7~T~T 

K=NCP 

□&= = — 

/ READ MODEL COORDS 

FOR ALL PTS ONE BY ONE 
















PRINT APPROXIMATE 
TRANSFORMATION PARAMETERS 


PRINT MODEL AND 
GROUND COORDINATES 



PRINT FULL GROUND 
VAR/COVAR MATRICES 


'PRINT LIST OF 
GROUND VARIANCES/ 



PRINT FULL MODEL 
VAR/COVAR MATRICES 


'PRINT LIST OF 
MODEL VARIANCES/ 


RETURN 
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4,4 Flow chart of INPUT3 



m 
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/ PR I NT VAR/COVAR MATRICES / 


/ FOR MODEL COORDS / 
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4.5 Flow chart of TRANSF 
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RETURN 
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THE FOLLOWING SUBROUTINES ARE INCLUDED TN THE PACKAGE AT PPFSEwT 

1, TNpuTl - FICTITIOUS DATA gENFRATOR for simulation 
?* Input? - reads d a ta from aerotri angulation outputs 
3, I NFUT3 - RFADS mOD f L COORDINATES AND v A R"CDvAR INFORMATION 

from aerctriangiilatipn OUTPUTS and simulates CONTRD 

LS FOR KFIATIVE ACCURACY TEST. 

A, TRANSF « LEAST SQUARE SOLUTION OF THE TRANSFORMATION 

5, ROTATE » COMPUTES POTATIDN MATRT y ANP ITS DErIvATTvFS WITH 

rFspect to the rotation angles 

6, GAUSS - COMPUTES NORMALLY DISTRIBUTED RANDOM NUMBr*S WITH 

GTvEn M E A n and deviation 

7* pandu = computes l'NIformely oIstributed random numbers 

8, weight - computes thf weights of the given points from the 

INpUT y A F*CO v AR m a TRTCFS AND ESTIMATED SIGMA 7FRD, 

9, DETCAL - COMPUTES thf CONTIBUTION OF each POINT TO I HE 

NORMAL EpliATlONS 

10, PART • COMPUTES THE CONTRIBUTION OF EACH POINT TU THF 

observation eouations 
n, mxmult - multiplies two matrices 

12, INVEPS " COMPUTES THF INVERSE OF A SouARE MATRIX BY DESTROY- 

ING THE ORIGINAL. 

13, sigma * COMPUTES The sigh* ZFRP Of the SOLUTION 
1A, RTODMS - TRANSFORMS radians to DFGREES-MIN-SEC. 

15, OUTPUT - AN OUTPUT SUBROUTINE FOR EACH ITERATION 
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C 16. rulPUF - FINAL output of the solution 

c 
c 
c 

c 

c 

c 

DIMENSION NUM( 100)#XOC 100)# YO( 100 >#Z 0 { tOO)#XOO( 100)# Y00( 100) 

1#700(100>#X( 1 00)# YUOO)* 7 noO)*WJ(iPO# 3# 3)* DON( 100# 3, 3 ) , DON I Ny( 100 
l,3,3),SlGj{100,3,3) # DDSlG(100,3,3),rDELTA(100,3),DOK(100,3),RN<100 
1,7#3)#SJGMJ(100#3#3)#DX(100)#DY(100)#07(100)#XMC100)# YM(100)#7M(10 
IO)#DfVlA<100#3#3) 

DIMENSION 0BT(7,3)#DDB1(3»3)#DPTWJ(7#3vDDPTMJ(3#3)#G(7,3)#DDk 1<3) 
l#0nK2(3)#P(7#7)#QQ(7)#BKi<7,3),BNT(3#7)#D0NIvT<7#7)#Hl<3#7)#H?( 3# 7 
1 >,H3(7#3)#A( 3# 3) , CB f 3/ 7 ) # DDB f 3 , 3 ) #C ( 3 ) # DDC { 3 ) ,DW ( 7 » 7 > # A W J ( 3 

|,3)#ADDw(3#3)#0DN(3,3>#rrNTNj(3#3)#rG(7)#Dn(7#7)#DC(7)#D r <K7),DDj( 
17#7),DM(7#7)#OK(7)#DOJNV(7#7)#DELTAC7),DHDC(7>#F<3#3)#SIG0ELC7#7), 
lDSlG(7#7)#C0R(7)#pAR(7)#SlGpARf7)#0FVPAR(7#7)#AS(3»3)#SlGPR(3)»SlG 
1M0(3)#V(3), jOEfi( 3 )#MJn(3)#SEC(3) 

DImFNSION AT{3#3)#MtJC0Upf 100# 3 ), GRCOOR f 1 00# 3 }# A TX ( 1 00# 3 )# TRANSL ( 3 ) 
C 
C 

C pASlC INPUT READING AND PRINTING 

0 

0 

RFAD(5#100 5 JJ#Np#NCP#MAylTR, jpRlNT, IpRlNT 
100 FORMAT ( 6110 ) 

RE AD(5# 1 1 0 ) SIGO 
110 FORMAT ( E ? 0 , 1 0 ) 

READ(5#150) ( SIGPAR(k)# K»l#7 ) 

ISO FDPMAT(7E10.3) 

DO 160 K s 1 # 7 
DO 160 L*l# 7 

JE ( K - L ) 161 # 16 ?, 16 J 
161 0SIG(K#L) s 0,0 
GO TO 160 

16? IE (K.EO.l) DSIG(K#L> ■ SIOPARfl) 

IF Ck»FQ#2) D5 1 G( K* L ) B SIGPAR(2) 

IF ( K ,F 0 * 3 ) DSIG(K#L) * SIGPAR< 3 ) 

IF (K,FO,«) DS I G ( K» L ) * SI GPAR (4 1 
IF CK.F0.5) DSIG(K#L) - S IGP AR f 5 ) 

IF (K,FQ, 6 ) 0SIG(K#L) = SIGPAR (6 J 
IF (K«F0»7) DSIG( K»L ) B SIGPAR(7) 

160 CONTINUE 
C 
C 

CALL OUTPUX <JJ#Np#NCP'M*xITR#JpRlNT,IpRlNT,SlGO#OSlG) 

C 

C 

IF(JJ»E0»1) GO to ?00 
IF(JJ.F0,2) GO to 300 
IF(JJ«F0»3) GO TO «00 
NRITE(6#500) 

500 FORMAT("!"#///////» 10*#" ERROR s rrONG PARAMETER GOVERNING lNpUT#( 
1 JJ)"#/# 10X#" CHECK THE FIPSt PARAMETER ON COLUMN 10 Of tHF FIRSt 0 
1 AT A CAPO") 

GO TO 1000 

200 CALL INpUT 1 fNCP, xO# YO » zO# X# Y# 7 , 5CAL# OME # PH I # C APA# XC# YC * ZC# NUM 

1,SIGJ#0DSIG ) 

GO TO 999 

300 CALL INPUT2 ( Np# N C P * xO, YD# ZO# X# Y# Z# S IGO, SC AL» Xl# YC# ZC# OME# P 

I HI* CAPA, jPRINT,IPRlNT,SIfj,DDSIG , NUM) 
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C 

C 

c 

c 

c 

c 

f 

c 

c 

r 

c 

c 

c 

c 

c 

c 

r 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r 

c 

c 

c 

r 

c 

c 


GO TO 999 

4 00 CAl L INHUT3(NCP,Niiy,X0, Y0, 70, Y , Y » Z, SC AL* OMl * PH I , C APA , XC, YC.7 

1C,0DSIG,SlGj) 

Gn I 0 999 

999 CALL TRANSr ( NC H ** P * NUM , XO, Y0 # ZO, X, V , l» SC AL» OMt , PH I * C Ap A , 51 

1GJ,DSIG,DDSIG>SIG0,JJ,MAXITR,JPRINT,IPPINT ,XC,YC,7C) 

1000 CONTINUE 
STOP 
ENO 

SUBROUTINE INPUT l<Np,xU,YU,z0,y>Y>Z»SCAL,0ME,PHI,CAPA»yC>fC,7C#NUM 
l,STGJ,PDSlG ) 


PURPOSE 

GENERATES model and GROUND COORDINATES Up gIv En STANDARD pE v i a t I dn 
fop A set Df IDEAL P0IN1S, THE TRANSFORMATION PARAMETERS A RF Al SO 

random! y generated and tpetr approximations are randomly chosen 


pescription or parameters 

1 . inputs NP - number tf points used 

GXO - 

GYO - LIST OF "TRUE" COORDINATES OF POINTS 
G 70 - 
SIGGRX * 

SlGGRY - VARIANCE Op GROUND COORDINATES OF POINTS 
SIGGRZ - 
SIGMOX • 

SIgmOy - VARIANCE Of MODEL COORDINATES OF POINTS 
SIGM 07 • 

?. OUTPUTS XO - 

YO - PERTURPED GROUND COORDINATES 

CONTINUE 

70 - 

X ■ 

Y - PfRTURpEO MODEL COORDINATES 

7 * 

SCAL - APPROXIMATE SrAlE Op THE MODEL 

xe - 

YC - APPROXIMATE TRANSLATIONS OF THE MODEL 

zc - 

OME - APPROXIMATE rotation ABOUT X Axis 

PHI - APPROXIMATE ROTATION ABOUT Y AxlS 

CAPA - APPROXIMATE rotation ABOUT Z AXIS 

remarks 

ThE SUBROUTINE WORKS 0 N l Y FOR UNIFORM GROUND AND MODFL ACCURACY 
IT ALLOWS FOR DIFFERENT ACCURACIES TN EACH OF THE THREE BASIC DI- 
MENSIONS FOR BOTH MODEL AND GROUND, 

SUBROUTINES REQUIRED 

GAUSS ( normally distributed NUMBER GENERATOR } 

RANDU f UNIFORMLY distributed numbep generator ) 

ROTATE ( COMPUTES ThE ROTATION MATRIX TOR 3 -D TRANSFORMATION ) 


i « « « • « « • I f I t A » I f | • M f M I I » f • f • « « • M • f I M M ft t 9 M * | t | « » f M M • I I 

01 MENS I ON GXO(lOO)*GYO(tOO)#GZO(lCO}»xn£lOO)»YO( 100 ),ZO(lOO)>x (100 
l),YLIOO)»z(100),XM(100i#YMfl00 5#7MMOOI#AC3*3),ASC3#3>#NUM<100)»DS 
lIG(7#7>,SlGj(100,3,3),DrSlG(!00,3,3),PAR(7),COR(7) P DX(lOO),DY(lOO) 

1*07(100) 


C 

C 


REA0(5*20) GSCALpIAREXP 
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20 FORMAT(FlO. 5*110) ■ 

NPITF(6»30) 

30 FORMAT ("l", w *********************** PRINT UUT UF THE INPUT **** 
*00 35 Ja 1 , NP 

35 RFAD<5*45> GX0{ J), GYP ( J ) » GZO< J ) 

45 FORMAT! 3F15.5 ) 

WRITF (6*50 ) 

50 FORRATf /////, 3x # f LIST Pf TRUE cOORplNATEs OF POINTS USED IN THE S 
10LUTI0N***/,l6X,** POINT*', 23X,*> COORD T NATES** ,/, 1 5 X #" NUKBf H", 1 3x,** X 
I", 12V#** T*», 1 2X, " 7") 

DO 60 Jal,NP 

60 WFITF(6,70) J*GXO( J)*GTDf Jl,GZO( J> 

70 F0RNAT(16x»I3» 1 2X, F9 . 3 * 2 ( 5X , F9 . 3 ) ) 

WPITE<6,75) 

75 FCRMAT<*»!«, ////,* ******************** VARIANCE-COVARIANCE MATRIC 
1 WRITF(6,80) 

80 FOR MAT ( 3X'** THE VAR I ANcE - CP V AR I A NcF MATRICES FOR MODEL COORDlN 

1 ATFS ***///) 

RFAD(5^85) SIGM0X,SIGMDY, SIGM07 
RFAD (5* 85 ) SIGGRX*SIGGKY*SIGGR7 

B5 FORMATt 3E 1 0 . 3 ) 

SGX a SQRT(SIGGRX) 

SGY a SQRT(SIGGRY) 

SG7 = St>PT(SIGGR7> 

SMX a SoRTfSIGMOX) 

SHY * SqRT(SIGMOY) 

SM 7 = SORT (SI GM07 ) 

RAT I Ox * SGX * ( GScAL/SMX ) 

RATinY a SC.Y * C GScAL^SMY ) 

RATI07 = SGZ * C GSCAL/SMZ ) 

DO 90 Ja 1 , NP 
DO 90 K = l,3 
DO 90 Lai, 3 
lf( K-l ) 91*92,91 

91 SIGjfJ*K*L) a 0.0 
GO TO 90 

9? IF(K,EC,1> SIGj(J,K,L) a SlGMOy 
IF( K «E©.2) S!GJ(J,K,L) a SIGMOY 
IF (K.E0.3) SIGJ(J,K,L) a SIGM07 

90 CONTINUE 

DP 95 Ja 1 , NP 

95 NUM{J) a J 
J * 1 

OP 100 K* 1 , 3 

100 WPITF(6,130)(SIGJ( J,K,L)#L=1,3) 

130 F0RMAT(8X,3(E10.3,5X)) 

WPITE(6,190) 

190 FCRMAT(///»3X#** THE V AR I A nC E*C 0 V AR I A NC F MATRICES FDR GROUND COORD! 
1NATFS **,///) 
on 200 J® 1 , N p 
OD 200 K*l,3 
DD 200 L* 1 , 3 
TF( K-L ) 201,202,201 

201 DPSIG(J,K,L) a 0.0 
GP TO 200 

202 IF(K.E0,1)DDS1G( J,K,L) a SlGGRy 
IF(K.EQ,Z)DDSIG{J,K,L) a STGGRY 
IF(K.E0,3)DDSIGCJ,K,L) a SIGGR7 

200 CONTINUE 

J a 1 
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OP 205 K e 1/3 

205 WR1TE{6,130) (DDSIG(J,K,t )»L»1,3) 

WRITE (f»2l0) RATIO* 

WRITE (6/211 ) RAHOV 
WPITEt/S/212) RATIO? 

210 FPRMATc//" THE RATIO $IpMA GROUND TO SIGMA MODEL / EOR T Ht X COOROI 
1 NATES/ AT THE GROUNO SCAIE ISi*,Fl5.7) 

211 FpRMAT C// " THE RATIO slGMA GROUND TP sTGMA MODEL/ EpR THE Y CflOROI 
1 NATE S/ A 1 THF GROUND SCAIE I5",F15.7> 

21? f-ORMATf//" THE RATIO $IfN* GROUND Tp STGM* MODEL/ fOR ThE 7 CpOROI 
1NATFS/ AT THE GROUND SCAIE IS",F15.7) 

C 

C 

c generate transformation parameters 

c 

c 

RFAD(5/?30)GXC/GYC»GZC/GDMF/GPHI /GCAPA 
230 FPRMATfl0X/6F10.5) 

C 

C 

c realign ground points 

c 

c 

READ(5/315t IX1/S1/IX2/S? 

315 EORMAT(?( I10/E15.5) ) 

AM s 0.0 
IX e m 
S a S 1 

DP 310 1= 1,NP 

CALL GAUSS (IX,S/AM,C0HCPX> 

GXDtn sCGXO(l) ♦ COFGHX )*10,**IAREyp 
CALL GAUSS (IX,S,AM.C0RGPY1 
310 GYO(I> r(GYOCl) ♦ C ORGR Y )* 1 0 , * * I ARE YR 
IX = 1X2 
S s S2 

DD 320 1= 1»NP 

CALL GAUSS ( IX/S/AM,CorGp2) 

320 G70(I) e(G70CI) + C0RGK7)*AO» 

C 

C 

C generate model coordinates 

C 

C 

CALL ROTATE f GOME/ GPHI / GC ApA, A/ AS T 
DD 340 I = 1 » N p 

ox( n = Gxncn-Gxc 
D Y C I ) e GYOC I J-GYC 
07 C I ) = GZPC I )“G?C 

XMU) « (Af 1 ,1 )*DXC| )*A(1 ,2)*0Y( I ) + A(l ,3)*DZ( I ))*GSCAL 
YM ( I ) * (AC2/n*DX<n + AC?/2)*DY(I) + A(?/3)*D7(n)*GSCAL 
3A0 ZM ( l ) e ( A( 3/1 )*DX( I )*A( 3»?)*0YC I) + a( 3/3>*nZ( I) )*GSC aL 


perturb the model COORDINATES 
c 
c 

RE ADC 5/ 345) I X3/ 1X4/ I X5/ I X*/ 1x7/ IXfl 
345 FORMAT C 6l 1 0 ) 

SMX = SqRtCSIGMOX) 

SMY s SQRTCSIGMOY) 

SM7 * SQRTcSIGMOZ) 

DO 350 I«1/NP 
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CALL 

GAUSS 

CIX3#SMX*AM*VMX) 

X C I ) 

* XM(I) 

♦ VMX 

CALL 

GAUSS 

(IX4#SMY*AM»VMY) 

Y( I ) 

* YM(I) 

+ VMY 

call 

GAUSS 

CIX5*SMZ*AM*VMZ) 

350 Z(I) 

= ZM( I ) 

♦ VMZ 

PEPTURP 

GRODNO 

coordinates 


DG 360 I*1*NP 
SOX s SQRT(SlGGRX) 

SGY = SOFT ( S I GGRY ) 

$G7 = SqHT(SIGGRZ) 

CALL GAUSS ( 1 X6* SGX* Ap, VOX ) 

XO(I) = GXO(I) + VGX 

CALI GAUSS (IX7/SGY>AK,VGY) 

YP { I } * GYO( I) + VGY 

CALL GAUSS (IXB*SGZ,AM,VGZ) 

360 ZO(I) e GZP(l) + VG J 
C 
c 

c choose approximations for the transformation parameters 

c 

c 

RE ad (5, 370 5 DSCAL,DXC>DYC*O7C»D0mE,DPHT#UCApA 
370 FPRMAT(7F 10.6) 

scal * GSCAL ♦ DScAL 
DrF s gome + dome 
PHT * GPHI + DPHl 
CAPA s GCAPA ♦ DCAPA 

XC s GXC ♦ DXC 

YC = GYc ♦ DYC 

ZC b GZC ♦ DZC 

WRITE (6*400) 

400 pORMATf *»1 "* " ********************* OUTPUT OF SUBROUTINE INpUTl * 
i ******************%//✓, iov** transformation par 
1 A M E T t R S"*/*10x>" - — 

1-- 5X*** PAPAMFTEr**,1?x*" GENERATFD"*5X*" CORRE C T I ONS", ?x 
1*" APPROXIMATIONS**# /) 

WRITE(6*410) GSCAL»DSCAl #SCAL 
410 F0RMATI5X#** SC ALE"> 17X* E 1 0. 5* 2( 5X *E1 0, 5 ) ) 

WRITE (6*420) GXC»DXC*XC 

420 FORMAT (5X* •* TRANSLATION IN X**# 5X> Ft 1 « 5* 2 <4X* Fl ! . 4 ) ) 

WRITE(6*430) GYC#DYC»YC 

430 F0RMATC5X*** TRANSLATION IN Y"*5X#FM ,5*2(4X*F1 1 *« ) ) 

HR I TE( 6* 440 ) GZC*DZC»ZC 

440 FORMA T (5 X# ** TRANSLATION IN Z** 5 X# F 11 » 5* 2 (4 X* F 1 1 • 4 ) ) 

NR I TF ( 6* 450 ) GOME#DOME*rME 
450 FDRMAT<5X#** AnGLF OME"' 1 3X » 3(F 1 0.5*«X ) ) 

WRITE (6# 460) GPHI * OPHI * PH I 
460 FORMAT! 5 X* * ANGLF PHI”# 1 3x» 3( F 1 0 , 5* 5 X >) 

WR I TE ( 6* 470 ) GCAPA* DC APA * CApA 
470 F0RMAT(5X*»* ANGLE C ApA"* 1 ?X * 3< Fl 0 , 5* 5X ) ) 

WR I TE( 6 * 475 ) 

475 FORMAT r/*5x*** ANGLES IN DEGREES'*) 

CALL RTODMS(C,OME#lDr.OME*MlGOME*SFGOME) 

CALL RTODMS(OOME*IL>EGDU*MIDOME*SEDOpE) 

CALL RTODmS{OmE*IDEgO*MINO»SECO) 
call rtoohscgphi^iogphi^migrh^segpmi > 

CALL RT0DMS(0PHI*IDEGDM*MIUPHI*SFDPHI) 
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C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r 

c 

c 

c 

c 


t 

476 


1 

477 


1 

4 7B 
4B0 

1 

1 

1 


490 

500 

5)0 

1 

1 

1 

5 ?0 


) 


call rtodmsc 
call rtoomsc 

CAU R T 0 U M S ( 
CALL RTOOMSC 
WRITE(6,476) 


PHI, IDEGP,MINP,SECP) 
GCAPA,IDGCPA>MlGCPA,SEGf P a ) 

DCAPA, IDFGUC,MTOCPA, SEDCPA 1 * 
CApA,IUEGC»MlNC,$ECC ) 

IDGPME,MIGOME,SEGOmE, If EGDU , M IDOME, SEDOM t , 


IUEG0,MINO, 


SEC 0 

FORMAT ( $X, " ANGLE OMEGA", 5X, 3 ( 2X , I 3, I 3* Y 6 , 2, IX ) > 

WRITE (6,477 ) I OGPHI , M I URN I , SFGPH I , I DEGDP, M I OPH I , SEDPH1 , 


SFCP 

FORMAT'S*," ANGLF PHI", 7X , 3{ ?X , I 3, I 3, F 6 . ?, 1 X ) ) 

WRITE (6, 478) ! DGC P A, Ml GC FA, SEGCP A , IT E GOC , M I DCPA , SEDC P A, 


IOEGP,MlNP, 

II)EGC»mINC, 


SECC 

FORMA T ( 5X, " ANGLF CAPA", 6 x , 3 ( ?X , I 3, I 3 , F 6 . 2, 1 X ) ) 
WPITE(6,480) 

F0RMatE////,?0X," POUF I. COORDINATE 

................ POINT", 14X, 

" PfRT»'RPeO",/,5X, m NUMPfR",7X," X*,9X," Y",9X," 

", 9 X , " 7 ",/) 

00 490 J«1,NP 

WRITE (6,500) J,XM(J),VM(J),2M(J),X(J),Y(JJ,Z<J) 
FORMAT' 7*, I3,5X,6f F9,2,?X )) 

WRITF (6,510) 

FflRMAT ( «! ",18x," GROl'Np COPRDlNAT 

................................ ,//, 6x ,h pm NT -,i 

3 *," PFRTURPE0",/,5X," NUMBER", 7 X," X",9X," Y", 9 X 


5",/,?0X," ----- 

" GENERATED", 23X, 
2", 9* , " * ", 9X , " Y 


E S-,/,lPx>" ---- 
4 X , " GFNFRATED",? 
," Z",9X," X" , 9X , 


" Y",9X," ?",/) 

DO 5?0 J=1,NP 

WRITE (6,500) J,GX0( J),GYO( J),G70( J),XO( J),YO( J),ZO( J) 
RFTliRN 


END 

SliRpnUTlNE INPUT? ( NR, FTP, xD, YO, 20, X, Y, Z, S I GO, SC A L, Xl , Yt , ZC, PME , p 
HI, CAPA, JPRINT, IPRINT#SIGJ»D0SIG ,NUM) 


PURPOSE 


1, REAPS MODEL COORDINATES LIE ALL POINTS (KNOWN AND UNKNUwN ) AT 
ANY SEQUENCE. 

?. reads the ground coordinates of the controls at any sEpuencE, 

3, READS VARI ANcE-COvAWI ANTE MATRICES ( FULL 3 y 3 MATRICES ) FOR 
Both GPoU n D and MoUEL C nr Rni NAT ES n F the CnNTRnLS AND nE (HE MqDFL 
COORDINATES of the UNKNOWN points. 

4, SEpArATFS CONTROLS ERDm UNKNOWN POINTS and puts The MODEL COORD 
INATFS IN The ORDER OF THF CONTROL POINTS. 

5, PUTS IN ORDER THE yApI AN cF-COvART ANcE MATRICES 

6, COMPUTES APPROXIMATION np THE TRANSFORMATION RARAmETFRS, 

DESCRIPTION OF PARAMETERS 


1. INPUT 


NP - 
NCP - 
M TO N A - 

NUM - 
* 0 , YO, zO " 
A, P ,C - 
SIGO - 
CONTINUE 


NUMpER OF model points 

NUMBER OF CONTNri POINTS 

numbers of the corner points of the model to ^e used for 
the computation of the transformation parameters. 

ID NUMBER OF THF POINTS 

GROUND COORDINATFS OF THE CPNTPULS 

MODEL COOROINAIFS Of TpE MODEL POINTS IN ARBITRARY ORpER 

estimated variance of unit weight 


SS1-SS9 - ElFmENTS OF THE VAR-COvAR MATRICES OF THE MODEL POINTS BY 
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C ROW IN ARBITRARY ORDER 

c si -$9 - elfments of the var/covap matrices of the ground coordIna- 
c tes nr the controls by row in arbitrary order 

C DSIG - VAR/COVAR MATRIX FOR THE APPROXIMATE TRANSFORM! ICN PARAM. 

2 , OUTPUT 

xn#Y 0 »zn - ground cdordinatfs for the controls 

x,y»? - model coords FOR all model points in correct order 
S lGJ - VAR/COVAR matrices of the model coords in correct order- 
c NP# 3/ 3 ) 

DOS I G - VAR/COVAR MatRKFS FOR jHE GROUND COORDS Of THE CONTROL 
POINTS IN CORRECT ORDER ( NCP, 3#3 ) 

scal - approximate scale of the model 

XC,YC, 7 C - APPROXIMATE TRANSLATIONS of thf mooel 
CONTINUE 

ome - approximate rotation uf the model arround the x axis 

pht - approximate rotation of the model arround the y axis 

capa - approximate rotation of the model arround the z axis 

REMARKS 

* ThE SCALE Is COMPUTED as the mean of the SCALES Of THE two 
DIAGONALS Of The model 

* THE TRANSLATIONS ApF COMPUTED AS THE DIFFERENCES IN x»Y> AND 
7 OF THE CENTFp UF GRAVITY OF An CONTROL POINTS 

* the four corner p uifts ni,n?,n 3 ,na arf used for the computa- 
tion of the rotations. 


DIMENSION YO(100)#YO(100v70(100)>XflOO)>Y(100),i!(100)>NUM(100)# D 
IDSlGf IDO' 3)P SIGJC 100' 3 p 3 ) 

read the ground coordinates of the control points 

READf5#10) N1*N2»N3,N4 
10 FORMAT(AllO) 

oo no i* i # nc p 

no RFADfSnOO) NUM ( I )»XO(I ) » Y 0 ( I )» ZD ( I 1 
100 F0RMATU5#3F15.3) 

read the model coordinates of the points 

K s NCP 

400 RF AD (5# 450 ) N#A,B>C,II 

A50 form AT CIS# FI 6, 3# FI A. 3#n 5, 3# 26 X* 1 2 > 

J e 1 

310 IF(N ,E0 ,NUM{ J) ) GO TO *C0 
IFCJ.GT.NCP) GO TO 300 
J s J + 1 
GO TO 310 
200 X(J) * A 
Y ( J ) » B 
ZCJ) * C 

350 IF(I!»NE«99) GO TO 400 
GO TO 500 
300 K « X 4 1 

NUM(K)«N 
X(K) = A 
Y(K) s B 
Z(K) * C 
GO TO 350 
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500 CONTINUE 

c 

c RF A[l AND PUT IN ORDER TrF VAR/COVAR MATKICFS OF 1 HF CONTROLS 

C 

600 FpRMAT(l5*9EB«2>l3) 

RE AO (5, 600} N>S1*S?»S3* S* > S5> $6, $7# s8, $9, 1 1 1 
J = t 

6?0 IFfN.EO.NUMtJ)) GO TO 6JP 
J s J + 1 

GP TO 620 

610 DPSIGCJ# 1>1 ) * SI 
0f)SlG(j>l>2) = S? 

DDS1C,(J,1,3) * S3 
DPSIG(J»?*n = SA 
DOSIGC Jj.2»2) * S5 
OPS JG( J,2>3) * S 6 
DDSIG(J,3,i) * S7 
DPSJG(J,3,2) * S 8 
DOS TG( J# 3# 3) = $9 
IF (TII.NE.99) Gfl TO 500 
C 

c REAP AN'd put IN ORDER IrF vAR/cOvAR MATRICES Of ThF model POINTS 

c 

K a NCP 

700 READ(5^600) N> $S 1 * S$2* S$3 * SS A* S55* SS6> SS7> SS8* SS9> II 11 
J * 1 

730 IFCN ,E<\NUMC JI) GD TO 710 
J = J + 1 
GO TO 730 

710 SIGJCJMM) a SSI 
SIGJ(J>1’2) * SS? 

SIGJfJM»3) * SS 3 
SIGJ(J#2*1) = SSA 
$IGJ(J*2*2) a SS5 
SIGJ(J#2#3) a SS6 
SIGJ(J*3*1) * SS7 
S!GjfJ,3*2) a SS8 
SJGJ(J>3»3) s SS9 
750 lF(ini,NE,99) GO TO 700 
GO TO 7A0 
7 A 0 CONTINUE 
C 

C COMPUTATION OF AppRDKl MATIfiNS FOR TuE TRANSFORMATION rArAmETFrS 

C 

C 

DO 900 J a 1 * NCP 
IFCNIJMC J),EQ,Nl) K a J 
IF(NUMCJ),F0,N3) L*J 
IFfNUMC J),F0,N2) Maj 
900 IFCNUMf JJ.FQ.NA) N*J 

computation of approximate scale 

DMl S AgS(SoRT{ (X (K )*X(l ) )**?+(Y(K)-Y(L ) )**2+(Z(K )-Z(L ) )**2) ) 

Dm2 * ABSCS0RT((X{MW(N))**2+(Y<M)-Y(N))**?«'{Z(M)-Z(N))**2>) 

DG1 rABS(S0PT((XO(K )-XO(l ) ) * *2 + f Y 0 ( K 5 -YP ( L ) >* *2 ♦{ ZO f K ) -ZD ( l ) ) ** 2 ) ) 
DG2aA8S(S0RT({X0fM)-X0(N5)**2+F YO(M)-YP(N))** 2 +(Z 0 (M J»Z 0 (N))** 2 )) 

SCAL « (fDMl/OGl) + ( DM2/DG2 ) )/2 
COMPUTATION OF APPROXIMATE TRANSLATIONS 



39 


SUMXM * 0, 

SUM YH e 0, 

SUMZM a 0, 

SUMXG * 0. 

SUMYG * 0, 

SUM7G * 0, 

DO 910 I * 1 * NCP 
SUMXM * SUMXM + X(I ) 

SUMYM = SUMYM + Y(I> 

SUMZM e SUMZM + 7(1) 

SUMXG s SUMXG + X0( 1 ) 

SUMYG « SUMYG + YOCI) 

910 SUMZG e SUMZG ♦ ZOCI) 

XC» (SUm>G - SUmXm)/NCP 
YCe(SUMYG - sumymj/mcp 
ZCe(SUMZG* SUMZM5/MCP 
C 

c computation of ROTATION CAPA 

c 

DXKL a X<K) - X(L) 

DYKL « Y(K) - Y ( L ) 

DXOKL = XO ( K ) - XO ( L 5 

DYOKL = YO(K) «■ YO(L ) 

DYOMN s YO(M) - YO(N) 

DXOMN s XO(M) - XO ( N ) 

OXMN e X ( M ) - X ( N ) 

DYMN a YCM) • Y ( N ) 

THETM1 a ATANZ(DYKL ,DXK» ) 

THFTG1 E ATAN2(DYDKL#DXDKL) 

ThFTmZ a ATAN2(DYMN*DXMN) 

TMFTG2 e ATAN2cDY0MN,DXrMN) 

CAPA s ( ThETgWThFTg2”TMFTh1-ThFTm2)/2 

c COMPUTATION OF ROTATION PHI 

C 

DzMK * ZfM)-Z(K) 

D7L N a ZCL)-Z(N) 

X YMK s S0RT((X(M)-X(K))**2 + ( Y ( M )-v ( K ) )* *2 ) 

X yI N b S0RT((X(L)-X(N>)**2 ♦ ( Y( L )»v( N) )**?) 

DZOMK a ZOCM)-ZO(K) 

D70LN a ZO(L)-ZO<N) 

XYOMK a SORT ( (XO(M)-XO(K ) ) **2 ♦ ( YO (M ) -YO ( K ) ) **2 ) 
XYOIN a SfiRT((XO(L)-XOtN3>**2 ♦ ( Y0( L ) - YOC N ) )**2 ) 
PHIM2 a ATAN2(DZLN*XYUN) 

PHI M 1 e ATAN2(DZMK# XYMK ) 

PHI G 1 e ATAN2CDZPMK, XYOMK) 

PHIG2 a ATAN2(DZ0LN,XY0t M) 

PHI = ( PHIG1 ♦ PHIG? - PHIM1 - PHIM2 ) / 2 
C 

c COMPUTATION of ROTATION omega 

c 

DzKN a ZOO • Z(N) 

D7ML a Z(M) - Z(L) 

XYKN a S«RT((X(K)-X(N))**2 ♦ ( Y ( K )-Y ( N ) )**2 ) 

XYMl a SGRT((X(M)-X(L))**2 ♦ ( Y(M)-Y(t) )**?) 
DZ0KN»70(K)-Z0CN) 

DZOMLaZO(M)- ZO(L ) 

XYOKN e SORT((XO(K)*XO(N3 )**2 ♦ ( YO ( K )-Y0 ( N ) )**2 ) 
XYOHL e SOPT((X0(m 3-X0(L))**2 4 (Y0(M)-Y0(l))**2) 
0MFM1 a ATAN?(DZKN,XYKN) 

OMEMP b AT AN2( DZNL* XYML ) 

0MEG1 e ATAN2(DZ0KN,XY0KK) 
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0MEG2 = ATAN2CDZ0ML.XY0MI ) 

OMF S ( UMEGl ♦ OMFG? -CMtMJ - OMEM? ) / 2 
CALL RTODMS(OME. lD0ME.MpNf.SOMF) 

CALL RTODMS(PHI.IDPHl.MPNl.SPHl) 

CALL RTODMS(CAPA, I OC A P A . NC Ap A . SC AP A ) 

C 

c print rui of the olitpui of subroutine input? 
r 

WPITE(6.1000) 

1000 FpRMAT( "l**.** ********************* OUTPUT OF SUBROUTINE INPUT? * 
i *******************. / / / . i OX . ** TRANSFORMATION par 

1A M F T F R S"#/.10x." * 

PARAMFTER",20x.** COMPUTED A PPRUx I MA 1 I DNS**, /, 3 gx . 

1" IN RAD**>9y»** IN DEGREES**./) 

WRITF(6#1010) SCaL.XC.YC.ZC 

join FORMAT ( 9 X. •* SCALE**.30X#no,5./. 9X.** TRANSLATION IN X **. X * X# F 1 1 . 3# 
1/.9X." TRANSLATION IN V " . 1 AX . F 11 » 3# /. 9X . ** TRANSLATION IN /".1BX.F1 
11.3) 

WRl TE ( 6 . 1 020 ) OPE. IDOMr.MOMF.SOME 
1020 FPPMATf 9X » '* ROTATION UpF G A " , 1 1 x ,F 1 0 ,5 . 9X .2 C 1 3 , 1 X ) * F 7 , 3 ) 
WPITE(6,1030) PH I . I DPH 1 . NPH I . S PHI 
1 030 FORMAT! 9 X.** ROjATIPn *• . 1 IX .F 1 r . 5. »X * ?( 1 3. 1 X ) . F 7 . 3 ) 

WRITE (6,1040 ) CAPA.IDCAFA.MCApA.SCAPA 
1 0 40 FOPMATf 9 X » ** ROTATION CAPA **. 1 1 X * F 1 0 . 5 * 9X. 2 ( I 3 . 1 X ) » F 7 . 3 ) 

WRITE (6.1100) 

1100 F0RMAT(////#20X»" MODFL COORDINATE S%/#20X»" 

1 ............................ **,//, 6X#** 1. CONTROL POINIS".//. 

16X." NUMB E R w » 5X . " POINT*. 10X." X".1CX." Y**. i OX . «Z*» » / ) 

DO 1120 Jsl.NCP 

1120 WRITF(6.1130) J.NUM( J).X( J).Y( J).Z( J) 

1130 FORMAT ( 5X.I5.8X. 15. 3F12. 4 ) 

IF(NP.FO.NCP) GO TP 1160 
WRITEC6.11A0) 

1140 FORMAT ( // . 6X . ** 2. UNKNUWN POINTS*’.//. 

16X.” Nt'MPf P*. 5X. W PPIN1 •’.lCX.* X".1CX.** Y OX ."/**./ ) 

NCPPl*NCP+i 
Dp 1150 JeNCPPl.NP 

1150 WRI TF (6, 1 1 30 ) J.NUM( J).X( J).Y( J),Z(J) 

1160 WPTTF(6.1170) 

1170 FPRMATCM". 18X." GROUND COORDINATE S**. /. lt*X» " — — 

6X.« CONTROL POINTS".//. 

16V.** N1MPEP".5X." POINT*. 10X." X**.10X." Y« . 1 OX .**Z" . / ) 

DO 1180 Jsl.NCP 

1180 WPITE( 6. 1 130) J.NUM(J).Xr(J). Ynt J).70( J) 

IF(JPRINT,F0,1 ) GO TO 1500 
WPI TF ( 6. 1 31 0 ) 

1310 F0RMAT(*1 *»3X.* THE VA RI A NCE "COVAR I ANfF MATRICES FOR GROUND COORDI 

1 NATES**.///) 

DO 1320 J sl.NCP 

1 320 WRITE(6, 1330) NUM ( J ) . ( ( DOS IG( J. K . I ) . K* 1 . 3 ) . L= 1 . 3 ) 

1 330 FORM ATf 3X.** PO IN T*M5. 2* » 3( 2X . E 1 3 , 6 ) . /* 2( 1 6X . 3 < 2X . F 1 3 • 6 ) . / ) ) 

GO TO 1340 
1500 WRITFC6, 1510) 

1510 FDRMAT(»1 ".///. 3y." THE l 1ST OF VARIANCES FOR GROUND COORDINATES". 

1//.5X." number**. 3x.** point*. 7x.** vap x**.8x.»* var y^bx.** var ?*•) 

DO 1520 Jcl.NCP 

1520 WPITE(6. 1530) J.NUM( j).(ri?SlG( J#K.K).K*1 .3) 

1 530 F0RMat(I10.5X#I5.5X.3(E)?.5,2X) ) 

1340 IF(IPRINT.EO.I) GO TO 17C0 
WRITE(6. 1610) 

1610 eormat(«i».3x.” the variance-covariance matrices eop model COORDIN 

1ATES **.///) 
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DO 1620 J sl,NP 

16?0 WRITE(6> 1330) NUM<J)*(( S I G J ( J , K , I ) , R* 1 * 3 ) , L» 1 # 3 ) 

GO TO 16«0 
1 TOO WRITE(6,1710> 

1710 rORMATf *\* >///>**>* tHL I I S T Of VARIANTS TOR MDOeL COORDINATES", 
1 //* 5X# ** KUMBER"#3X#" POINT", 7X," VAP X"»8X*" VAR Y",RX," VAR 7") 

DO 1720 J* 1 , NP 

1 7?0 WRITf(6,1530)J>NUM{J),( SIGj(j,K»K),Kb1,3) 

1 6 AO CONTINUE 
RETURN 
END 

SURROUTINE INPUT3<NCP,NUP>X0# YO*ZO* yp Y#Z#SCAL#OME>PHI#CAPA,xC, yc>z 
1C,00SIG,SIGJ) 

c 

c 


this subroutine reads uuTput from any triangulation program c mo- 
ofi. coordinates AND THEIR DEVIATIONS )0 

it rotates the model by small rotation angles and the computed 

COORDINATES ARE cONslD^D TO pE FICTITIOUS "GROUNO" CONTROLS, 

the c, round is assigned arbitrarily nigh weights and the program 

COMPUTES The ADJUSTED TRANSFORMATION PARAMETERS ano THEIR DEVIA- 
TIONS, WHICH ARE REPRESENTING ThE UNCERTAINTY OF ThE pOOEt AS A 

whole, 

ThE INPUT COORDINATES Art DEVIATIONS need not bE IN TrE SAmE order 


DIMENSION NUM(100),X0(100>, YOf lOO),ZOf 100 i, xc 100 Yf 100 W(iOO), 

1 D0STGfl 0 0,3'3),$lGj(10C,3j , 3),TDEG(3),MlN{3),SFc(3)>A(3#3),AS(3>3) 

1 , V( 3) # S 1GGR( 3) 

READ INpUT 

READ(5>5) IX1>IX2/IX3^IXA 

5 FORMAT(AIIO) 

DO 10 J=1#nCP 

10 READ(5,20) NUM(J), XCJ), Y(J), Z(J) 

20 EORMATf I10,3E15,2) 

35 I a 1 

REA0f5,40) N# S 1 , S2* S3* 1 1 1 1 T 
40 FCRMATC llO*3F20, 1 0#6X» 12) 

45 I E ( N , EO » NUM { I ) ) GO TO 30 
I e I + 1 
GO TO 45 

30 $IGJ(I,i*l) ■ SI 
SIGJ(I>1#2) a 0,0 
SI Gj( 1,1*3) = 0.0 
SIGJCI*2#1) * 0,0 
SIGJ(I*2»2) b $2 
SIC,JCI»2#3) b 0,0 
S I G J( I * 3# 1 ) b 0,0 
SIGJC I »3»2) e 0,0 
SIG J( I* 3* 3) « S3 
IE(IIIII.NE,99) GO TO 3& 

PRINT THE READ INPUT 

WPITEC6,50) 

50 fORMAT(" 1",///,5X>** THE PlvEN MODEL COORDINATES OF ThE CONTROL PO 
i INTS"*///*SX»" NUMBER"*! 1 X",!TX,** Y”#17X#** Z"*/) 
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00 60 T = 1 ,NCP 

60 WPITF(6,70) I,NUM(I), *(I), T(I), Z(I) 

70 F0RMAT(2l5,9x*3(FU .2,6y) ) 

WRITF(6,60) 

60 FORMAT! "1 ",///, 5y," THE OFvIATTONS PF THE MODEL CnORL-lNAIFS OF Th 
IF CONSOL POlN T S",///, 5X," NUMBrR",8X," SiGMa X",9X," SIGMA Y",9x 
1," SIG^A Z") 

DP 90 I = 1 , NCP 

90 WrITF( 6,100) I*NUMU),< ? ! G J( I , K , K ) , Kel * 3 ) 

100 FORMAT (2l5,9X,3(Fl 1 , 7, 6X ) ) 

C 

C GENERATE ARBITRARY ROTATIONS 

0 

RFAn(5,l05) DEG1,DEG2 
105 FORMAT ( 2F 10*5) 

IX s 1X1 

AM * DFG1/57. 2957795131 
S * OF 02/57*2957795131 
on 110 1=1,3 

CALL GAUSS! IX,S,AM, VC 1 )) 

110 CALL RTOPmSC v(I)*IOEG(I)»MlNCn,SECf 1)1 
OMF = V(l) 

PHI s V ( 2 ) 

CAPA = V f 3 > 

WPITF(6, 1 20) 

120 FORMAT!")", "*********************** PRINT OUT 01 INPUT3 ******** 
1 ***************„, /////, 5X," ThF ARBITRARY POTATIONS 01 THF MODEL", 
1 // , 30X , " IN RA[)«,9X," IN DEgREES",/) 

WPITF (6*130) V(n>IDEGfU,MIN(l), FtCfl) 

130 FORM A T ( 9X * * OMFG A", 1 3x* F 1 0 , 7, 7X* ?( 1 3* 1 X ) * F 7 . 3 ) 

WPITF (6,131) V(2), IDEG(2),MIN(?),SEC(2) 

131 FORMAT (9X," PHI ", 1 3X, F 1 0 . 7, 7 X, 2 ( P, 1 X ) , F7 . 3 ) 

WPITF (6*132) V(3),IDEG<3),MIN<3),$EC(3) 

132 F0RMAT(9X," CAPA ", 1 3 x* F 1 0 , 7, 7 X, 2 ( I 3 , 1 x ) , F 7 . 3 ) 

C 

CALL ROTATE (OME,PHl, CAPA »A, AS) 

C 

DO 170 1=1, NCP 

XO(I)= A( 1,1)* X(I) ♦ AM,?)* Y( I ) + A( 1,3)* Z(I) 

Yorn= AC2,n* xrn + Ar?,?j* ven * ac 2,3)* zrn 

zncn = a(3,d* xm + ao,?)* ycd + a( 3,3)* zri) 

170 CONTINUE 

c 

c asIgn values to approximate transformation parameters 
c 

RFAD(5,1C5) ATRA,DTRA 

IX =1X2 
AM S ATRA 
S = 0TPA 
DO 200 1=1,3 

200 CALL GAUSS(IX,S, AM, vUn 
XC = V( 1 ) 

YC = V ( 2 ) 

ZC = V ( 3 ) 

IX e 1X3 

READ(5,1UJ>) DEG3 , DF G4 
AM = DFG3/57, 2957795131 
S = OFGA/57 *29577951 31 
DO 210 1*1,3 

CALL GAUSS( IX,S,AM,V( I) ) 

V(I ) = - V C I ) 

210 CALL RTODMS < V (I ) , I DF G 1 1) , Ml N( I ) * SEC ( I) ) 
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o^r = V(D 
PHI » V(2> 

CAPA * V { 3 ) 

IX = 1X4 

RFAD(5, 105) ASC AL,USCAL 
AM r ASCAL 
S e OSCAl 

CALL GAU&S ( IX*S,Am,SCAL 1 
WR I TF ( 6, 220 } 

220 FORMATf ////////, 5X,« APPROXIMATION* FOR jHF TRANSFORMATION PaRAME 
1 TERS** //* 10x* * PARAmETFR ,, *24x* * APPROX ! M A T I ON", /* 38y, * IN RAD** 9x 
1*” IN DEGREES***/) 

WPITF(6,230) SCAl,XC, YC,7C 

230 FORMATC 9X,* SC ALE", 30X, F 1 0 , 5, /, 9X,** TRANSLATION IN X"# H*X»Fl 1 .5, 
1/,9x*" TRANSLATION IN Y*, 18x*fll .5* /* 9X*" translation IN /**1RX*F1 
n ,5) 

WRTTF(6,231 ) V(l)» IPEGCn*MlN<U#$FCm 

231 fORMATf 9X," ROTATION UpF G A"* H X , f 1 0 , 5* 9X , 2 < I 3, 1 X ) , f7 . 3 ) 

WRITER, 23?) V(2),inEGf2),MIN(2),SEC(2) 

23? FORMATS 9x*" ROTATION PHI "* 1 1 x* F 10 .5, 9x* 2U 3* 1 X > * FT , 3 ) 
WRITF<6,233) V(3>* irEG(3)*MlN( 3)*SfC<3> 

233 FORMATC 9X," ROTATION CAPA "* 1 1 X * F 1 0 , 5» 9X# 2 ( 1 3, 1 X )* FT * 3 ) 
WPITF<6,160) 

IPO fORmAT(«1*,////,5x*** GKCl’Nn COORDlNCTFS - (GENERATED bY ROTATION 0 
if the Given model coordinates)*, ///,sx*** mumrfR"»i 

11*,* X",17X," Y"*17X*" 7 "*/) 

00 190 I * 1 # NC P 

190 WRITE (6, TO ) !,NUM( 1 )* X0( I ), Y0( I ),ZOC I ) 

READ (5*250) < S! GGR < K) * 1 * 3 > 

250 FORMATf 3<E10,3*5X)> 

00 260 jsl , NCR 
DO 260 k s 1#3 
DO 260 L* 1 * 3 
IF (K-l) 261,262,261 

261 D0$IG<j*K*L> * 0,0 
GO TO 260 

262 DDSIG(J,K,L) b SIGGR(K)**? 

260 CONTINUE 

DO 160 Jel,NCP 
DO 160 K e 1 , 3 
DO 160 L e l,3 
IF (K-t ) 166*167*166 
166 SIGJ(J,K,U * 0,0 
GO TO 160 

16T SJGJ(J,K»L) = SIGjf J,K*K)**2 
160 CONTINUE 

WRITF(6, TOO ) 

TOO FORMAT f "1 *, ///," ******************** variance-covariance matric 
WPITF(6,T10) 

T10 FORMATf 3x*" ThF vA R I A NcE -CO VARl A NcF MATRICES FOR m 0D*-L COOROlN 
1ATES***///) 

DO T20 J B 1 , NCP 

T 20 WpITE(6,T30) NUMf J ), C C SlG J( J, K , l) * Kb 1 , 3 ) , l* 1 , 3 ) 

T 30 FORMATf 3X,** POI NT** 1 5* 2X* 3 f 2X , E 1 3 . 6 ) , /, 2( 1 6X* 3( 2X * El 3 , 6 ) , / ) ) 
WRITEC6,T50) 

T50 FORMATf ///*3X,* THE V A* ! A NO E-CO V AR I A NC F MATRICES FOR GROUnD COOROl 
1NATFS***///) 

WRITE(6,T60) (( P0StG(J,K,1),Kbi*3),L*1*3; 

T60 FORMATf 3( 1 6X* 3( 2X, E 1 3 . 6 ) , / ) ) 

RFTURN 

FND 



c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 
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c 

f 

c 

c 
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SUPROUTINE TRANSF C NCP#NP#NUM#XO#Yn#70#X#Y#Z#SCAL#OMt#PHl#CApA#SI 
1 Gj#DSIO#UDSlG#SIGO# JJ#MA*ITR, JPR I NT, I PR I NT #XC# YC#ZC ) 


PURPOSF 

PERFORMS ThF LEAST Spl'ARE f t t of onf digital model to AN other 

for REASONS OF CONVERSION ThE onf MODEL IS CALLED "MODEL" AND 
THE OTHER "GROUND", 

DFFINITIUN OF VARIABLES 

N p# N C P - I r T A L NUM^FR OF PUINTS AND NijmrLr Of CONTROLS 
HUM - 10 NIIMPER OF POINTS 
XD# YO # 70 - fi POUND COORDINATES of The CONTROLS 
X#Y#Z - MOTEL COORDINATES OF ALL POINTS 

scal - iff approximate scale of the model 
xc#yc#zc - approximate translations of thf model 
ome,pni,capa • approximate rotations of the modfl 

sigo - estimated mean square error of unit weight 

SlGj(Np,3,3) - VAP/CUvAR mATpIcfS Of The MODEL COORDS 
UF ALL THE POINTS 

DOSl G(NCP» 3, 3 ) - VAP/COVAP mATpIcFS Of THF GROUND COORDS 

of The controi points 

DSlG(7#7 ) - VAP/COvAR MATpIx Up T hE T RAN$f OR MA TI ON P*RA- 
MfTFRS 

MAxITR - MAXIMUM NUMpEP OF ITFRATIQNS AlLOwFP 


DIMFNSTON NUM(iOO),XO(lOO) # YOfl 00 ) # Z0( 1 00 ) # X00( 1 00 5 t YOOf 1 00 ) 

1 #Z00( 100>'*( 1 00)> Y M00>*?E1 0O)#WJ( 100# 3* 3) #DDW( loo* 3# 3)#D0NINVC 100 
1 »3#3)#SIGj(lOO#3#3)#ODSir(lOO#3#3)#PDFI TA( 1 00#3)#DOK(100, B)#RN(100 
1 #7#3 5 #SIGmjC 100#3#3),DX( 1 00)#OY( 100l#07( 100}#XM{ 100)# YMC 100)#7M{ 10 
10)#DFVIA(100#3#3) 

DIMENSION DBT(7# 3)#D0BI( 3a3)#DRTWJ(7# 31*D0PTW,J( 3# 3)# 0(7#3)# DDkK 3) 
l*DDK?C3)*R(7#7)#0O(7)#BNJ(7#3j#BNTE3*7)*00NlVTf7*7)#Hl(3*7)#H?C3#7 
1)»H3(7#3)#A(3#3> # 0 P < 3# 7 ) # ODB ( 3# 3 )# C ( 3 ) # DDC ( 3 ) . DW ( 7 # 7 ) # a W J ( 3 

l#3)#A00w(3#3)#DDN(3#3)*rrNlNj(3#3)#rGCn#D0(7#7)#0C(7)#0GJ(7)#D0j{ 
1 7#7)#DN(7#7)#DK(7)#0D!NV( 7#7)#0ELTAf 7)#DN0C(7 J #F ( 3 # 3 ) # S I GPEL ( 7# 7 )# 
lDS!GC7#7)#C0R(7)#pAR(7)#SlGPARf n#0FVpARC7#7)*ASC3#3)#SIGGR(3)#SlG 
1MO(3)*V(3),IOEG(3)#MINC3)#SEC(3) 


ASlGN APPROXIMATIONS TU THE TRANSFORMATION PARAmETFRS FOR THE 
FIRST ITERATION 

SCALO e SCAL 
XCO * XC 
YCO s VC 
ZCO s ZC 
OMEO = OME 
PHIO = PHI 
CAPAO = CAPA 
DO 5 J»1#NCP 
XOOCJ) e XO(J) 

YOD(J) b YO(J) 

5 ZOO(J) s ZO(J) 


C 

C 


SIGOLD » SIGO 
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c compute the weight matrices 

c 

c 

CALL WEIGHT < NP# NCP* J J»NUM* JPP I NT/ IPPI NT, S IGO, SI G J> 0SI G, DOSI 

1G,SCAL,WJ»DW,DDW ) 

C 

C 

ITEPsl 

C 

C 

C INITIALIZE matrices od AND pg 

c 

c 

470 OD 10 K S 1 , 7 
DG(K5 * 0,0 
DP 10 L*l,7 

10 ODCK»L) * 0.0 
C 
C 

c computation or elements Of rc ( right hand for transformation parAme- 
C TFrS OBSERVATION equations ) 
c 
c 

DC C 1 ) = SCAL-SCALD 
DC C ? 3 * XC-XCO 
DC C 3 > * TC-YCO 
DC C ft 3 = ZC-ZCO 
DC C 5 3 ■ OMF-OME0 
DC C 6 3 = PHI-PHIO 
DC ( 7 3 = C AP A**C AP A 0 
C 


compute The ROTATION mATr! y ANO ThE DEPlvATlyE FOR PHI 


CALL ROTATE ( OmE, PHI > C A pA . A, AS ) 


COMPUTE The CONTRIBUTION Of EACH POINT To The do and DG MATRICES 

c 

c 

DO 20 J=1 , NCP 

c 

c 

CALL PArT( A,A$,xO(J),YP(JWP(J)>XOO(J)>YOO(J),ZnO<J),XC,YC,ZC, 
lSCAl ,OP,DOB,C,DDC,X( J), Yf J),Z( J) ) 

C 

c 

00 30 K*1 , 3 

00 30 L* 1 > 3 
AWJ(K»L ) *WJ( J,K,L) 

30 ADDW(K,L> * DOWCJ,K,L> 

C 

C 

CALL OFTCAL (NP, j,DB,DOB,PC»ODC,C,AWj,ON,ADPN»ON,DDN,DK#DD«,BN, 
1D0NINV#0PNINj) 

C 

C 

DO 40 K* 1 , 3 
DO 40 L*1 , 7 

40 BNT ( K,l ) * BN( J,L»K) 

DO 50 K* 1 # 7 
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on 50 L = 1 » 3 
50 Bm J(K,I ) e Bw( J>K,l ) 

C 

CALL MXMULT(BNJ#DDNINJ#C*7»3#3) 

CALL WXMI»LTC0/BM,R,7# 7,3> 

DO AO K*1 »7 
OP 60 L B 1 # 7 

60 DPJ(K,l) c DN(K,L)-P<K>L ) 

00 70 K«l,7 

0 0 ( K ) s 0 # 0 
DO 70 L B 1 > 3 

70 OQ ( K ) * ©0 < K 3 ♦ 0<K#L)*PrKf J#L) 

DO 80 K* 1 * 7 

80 DGJ(K) = OKfK) «• QP(K) 

C 

c 

C CONTRIBUTION OF POINT J TU normal EQUATIONS 

C 

c 

DP 90 K*l,7 

DG(K)BrG<K)+OGj(K) 

DO 90 L e l>7 

90 DPCK>L)sPD(K,L)+DDJ(K,L) 

C 

?0 CONTINUE" 

C 

C 

c formulation of normal Equations 

c 

c 

DO 300 K=i,7 
DWOCCIO 8 0,0 
DO 300 LM>7 

300 DKDC(K) e DWDC(K5 + ON ( K# L T *PC (L ) 

DP 310 k*1»7 

DGCK ) = PG(K ) + PWDC (K ) 

DO 310 LM*7 

310 DD(K»L) b PD ( K » L ) ♦ DN ( K# L T 

C 

C 

c solution oe the normal Equations to obtain corrections eor the 
C TpANSfPRMATIpn pArAmETFRS 

C 

c 

DO 320 k* 1>7 
DO 320 L«l»7 
320 DOINV(K,l ) rPD ( K > L ) 

C 

C 

CALL INVARS ( DD I N V# 7# 1 1 ) 

C 

C 

DO 330 K*l*7 
DtLTACK)=0.0 
DO 330 L*l ,7 

330 DFLTA(K) e OELTA(K)+DDlNV(K>L )*DG(L) 


COMPUTATION Of CORRFcTIONS for thf ground coordinates op EACH point 


do 340 J*1 # NCP 
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Dfl 340 K«W3 

SUM 1 *0 . 0 
SUM2eO.O 
00 350 H*l*3 

SUM!* SUMl+DDNINVf J*K#M)*D0K( J*M) 

00 350 N*l»7 

350 sUM2»$UM2*DDNINV( J*K#M)*Pn( J» N# W ) *DFLTA < N > 

340 DDELTAC J# K ) * SUM 1 -SUM2 
C 
C 

C EppOR ANALYSIS 

C 

C 

CALL SIGMA (NCP*XO» YC#ZP>XOO* YOO*ZOP*xC'VC#ZC*SCAL#Df LTA jPOFLTA* 
1 A#AS>X#Y>Z#De»W»J#DDx#DN»SlGNFW ) 

C 

C 

DO 360 K B 1 # 7 
00 360 l*\»7 

360 S I GDELf K»L ) * S I GNE W+DUIN V ( K> L ) 

C 

c computation op tmF v A « IAN cF-covARiANcF matrices FOR EACH POINT 

c 

DO 410 Jb 1 # NCP 
DP 370 K e 1 » 7 
DO 370 L-l#3 

370 BNT ( L* K j = BN(J#K*L) 

DO 380 K e l>3 
DO 380 M = l#3 

380 DONI N J( M# K ) c DDN JNv ( J» M# K > 


CALL MXMULT < DDNIN j* BNT # HI , J* 7, 3 j 
CALL MXMULT { Hl#DDINy>H?#3*7*7 5 
C 
C 

DO 390 K B 1*7 
DO 390 L*1 » 3 
390 H3{ K # L ) = H1(L*K) 

C 

C 

CALL MXMULT { H2* H3» F # 3# 3# 7 ) 

C 

C 

DO 400 « B 1 * 3 
Dp 400 L B 1 * 3 

400 S I GM J( J# K # L ) c SI6NEw*(D0NlNJ(K#Lj<‘F(K#L)) 

410 CONTINUE 

ASI6N * SCRTCSIGNEW) 

Do 405 K B 1 # 7 

405 DEVPAR(K'K) * SGRTCSIGDEl ( k,k )) 

CALL RT0DMS(0EVPAR(5^5)» !DOEvO#MlDEvO#SEDEvO) 

CALL RToDMS(DEVPAR(6#6),!DOEVP,M!DEVPySEDEVP) 

CALL RT0DMS(DEVPAR(7» 7 1* I DOE VC# Ml DEVC* SEDE VC > 

C 

C 

C CONVERGENCY CRITERION 
C 

C IF(ABS(SlG0LD-SIGNEWl-SlG0*l,E-2> 420*430# 430 

430 CALL OUTpuT { I UR , ASI GN# DFLT A * Sf AL# XC* YC* ZC * OmL * PHI * C ApA, 

1DEVPAR ) 



ITFR * I T E R ♦ 1 

IF ( ITER-MXITR ) 440/440.450 
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C 

C 

c update VARIABLES 

c 

c 

MO Sf At =SrAL + PELTA( n 
XC sXC +DELTA(?) 
YC = YC +0ELTA<3) 
7 C sZf +DEL T A ( A ) 


OMF *OME +OELT A(5 ) 

phi sphi ♦pelt a ( 6 ) 

CAPA*CApA + pELTA(r ) 

DO 460 JM.NCP 
Xp(J) e XOCJ) + ODFLTACJ.n 
YOCJ) = YOfj) + DDE LTA ( J. ? ) 
460 ZPCJ) * Z Of J ) + DDE|TA(J.3l 
SIf, OLD « SlGNEw 
GO TO 470 


450 WPITF(6>640) HAXITR 

640 FOPMATf"!". 4 (/////// J / ?C X . ** 


1 /** * 
1?0X. 
t " 


>?7Xp" *"./.20x* * SOLUTION PDFS WOT CONVERGE *»./. 
* j AFTER". 15." ITERATIONS *"#//20 X/" *".?7X/ 


/ // 20 X 


GO TP 1000 

4?0 CALL OUT PUF f I Tt P. A S I GN» OF L T A . SC AL. XC / YC# ZC * DmE / PH I . c A pA. PF 

1 VpAp.SYGDEL.XO, YO.ZO.DDFL TA.NUM, NCp.STGMJ) 

IFfNP.F0.NCP5 GO TO 1000 

E ALL PTCOMP c Np.NCP.SCAL/XC.YC.7C/QME/PHI.CApA.X.Y/Z'NnM. 

5 PEVPAP.SIGJ) 

1000 CONTINUE 


RFTURN 


C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


END 

SUBROUTINE WEIGHT ( NP, NCP / J J/NUM# JPP INT * I PRI NT , S I GO# S 1 G J# DSI G. DDS I 
1G.SCAL.WJ/DWtDDW ) 


oo £ >aeooooooC>DOoQOOQDQOeooGeggQ I ?oc©ee&a** #t *e 9 aafied»*« # t t# *«*fl* ####t#t 

This SUBROUTINE COMPUTES TpE wFIfiHT MATRICES FOR ThE COORDINATE $ Or 
All TrF GROUND CONTROL POINTS Anp The MDpEL - CONTROL AND UNKNOWN - 
POINTS AND FDR The APPROXIMATE TRANSFORMATION PARAMETERS,. IT ACCEPTS 

full vari ancf-covariance matrices. 


INPUTS 

Np. NCP 
NUM 

JJ# JPPINT.IPRJNT 
SCAL 

Si go 

SlGj(NP.3,3) 


DDSlGfNCP. 3/3) 


DSlGfZ/Z ) 


iotal number of points and number of control 
id NUMBER of points 

INPUT PAPAMETFRS { SFE MAIN 5 
APPROXIMATE SCALF of the model 
estimated mean souarf error of lNIT wf i ght 

VAP/COVAR MATpICFS OF THE MODEL C00PDS 
OF AIL THE POINTS 

VAP/COVAR MATRICES OF THE GROUND COORDS 
OF THE control POINTS 

vap/covar matrix of the transformation para- 
meters 


OUTPUTS 


CONTINUE 

WJ(NP/3/3} - The WEIGHT MATRICES OF THE MODFL CUDROS OF 
AL I THE POINTS 
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C DW(7#7) - THF WEIGHT MATRIX Of THE TRANSFORMATION PARA 

MfTERS 

0DW(NCP»3#3) * IMF HEIGHT MATRICES OF THE GROUND tfUROS OF 
THE CONTROL POINTS 


OIMFNSION DW(7#7)#DSIG(7#7>#SIGJIN(3#3)*DDSINV{3#3)#dSIG 1P<7#7) 
DIMENSION Wj(lOO»3#3)»Up>‘nOO,3,3),rD$IG< 100#3#3)#SlGo< 10P#3#3)#NU 
1 M ( 1 00 ) 

C 

C INITIALIZE the WEIGHT MATPICE5 
C 

DO 100 K*l#7 
DO 100 L*1 #7 
100 DW( K #L )*0 « 

DO 110 Jel#NP 
DO 110 K e l # 3 
DO 110 L*l>3 
110 WJ(J#K*L>=0, 

DO 115 J®1 » NCP 
DO 115 K* 1 # 3 
DO 115 ltl,3 
115 DDW(J#K#L)*0. 

C 

C COMPUTE WEIGHTS or ThF TRANSFORMATION parameters 

c 

DO 200 K®1 # 7 

200 DWC*#K> * SIG0*<1/DSIG(K*K)) 
c PRINT wATrI y f OR TRANSFORMATION PARAMETERS 

c 

WR I TE { A# 1 45 ) 

145 FORM A T^"l ,, #////* ,, ********************** WEIGHTS USED IN SOLUTION 

220 f0RM^T(////' 3 X»" The weight MATRIX rOR The TRANSFORMATION pARAMETE 
IRS ”»///) 

DO 230 K * 1 # 7 

230 WF I TE ( 6# 240 ) { D W ( K> L )# L *1 # 7 ) 

240 F0RMAT(2X>7E10.3) 

C 

c COMPUTE weights for GPOUFD and model points 

c 

IF (JPPINT,NE,0 ) GO TU 400 
DO 320 J*1#NCP 
DO 310 K® 1# 3 
00 310 L e l * 3 

310 DOSINVHOL) b DDSIG(j,K,l ) 

CALL INVERS (DDSINV,3#II) 

DO 320 K*l>3 
DO 320 L«l*3 

320 DOW(J*W#l) b SIG0*DDSINVCK,L) 

GO TO 500 

400 DO 420 J®1*NCP 
DO 420 K®1 » 3 

420 DDW(J#K#K) e S I GO* ( 1 /DOS I G ( J#K* K ) ) 

500 IF CJJ.E«U2 ) GO TO 550 
J S 1 

WRITE (6*520 ) 

520 FORMATT////# 3X# " THE GpTUND CDORDlNATF WEIGHT MATRIX”#///) 

00 525 K* 1 # 3 
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5 25 WRITF(6,530) ( DOW ( J, K , L > , L * 1 , 3 ) 

530 F0RMAT(8X, 3( Ell ,«,5y) ) 

GO TP 56? 

550 WPITF(6,555) 

55S F0RMAT(////*3X,« THF WEIGHT MATRICES FOR THE GROUND COORDINATES Of 
1 THE CONTROL POINTS",///) 

J = 1 

5 70 WPTTF{6,560) NUMf J ) , ( DU* f J, 1 , L ) , l b 1 , 3 > , ( DDW ( J, 2, L ) , L a 1 , 3 ) , ( DD W( J, 3 
1,1 ),l =1,3) 

560 FnRNATf3X^ w POI NT"M 5, 3( F U , 4 * ?X ) * /, 20 AX, 3 ( E 1 1 . A* ?X ) * / ) ) 

IF(J-NCP) 561,562,56? 

561 J = J + 1 
Gp Tn 570 

562 iFf IPRINT ,NE*0) GO TO BOO 
00 720 j = J , NP 

DO 710 K B 1 , 3 
DO 710 1*1,3 

710 SIGjTNfK'L) * SIGj(j,K»LT 
CALL INVERS (SIGJIN,3,II) 

DO 720 K=l,3 
00 720 L = 1 , 3 

720 W J ( J, K , L ) = GlGo*sIGJIhfK,L) 

GO TO 900 

800 DO 820 j * 1 , NP 
on 820 K B 1 , 3 

B20 WJ{J,K,K) s SIG0*(1/SIGJ(J,K,K)> 

900 IE(JJ,NE,1) GO TO 1100 
1000 WR I TF ( 6, 1 020 ) 

1 020 FORPATc ////, 3 X, * TpF MODEL COORDINATE WEIGHT MATRIX"*///) 

J * t 

DO 1025 K»l,3 

1025 WRITFC6,530) C W Jf J* K* L > * L = 1 * 3 ) 

GO TO 1130 
1100 WRITF(6,1I?0) 

1120 pORpATf ///*3X*" THE w^GPT mATRIcES fOR THE COORDINATES Of THF MOD 
1 DEL POINTS f CONTROL AND UNKNOWN) ",//// ) 

J e 1 

IHO WrITF(6,560) NUMC J), f W Jf J, l,L ),L*1* 3 ), ( WJ( J,2,L)*L = 1 , 3), (WJC J, 3,L ) 

1,1 e 1 , 3 ) 

IFfJ-NP) 1125,1130,1130 
1125 J=J*1 

GO Tn 1 1 AO 
1130 CONTINUE 

rfturn 

END 

SURPOUTINE PART( A,AS»XO, Y0,zn,x00, Y0n,Z0O,XC,YC,7C,SCAL,0R,0DB,C 
1 , ODC, X* Y» 2) 

DTMFNSIqN DB(3,7),0DP(3,3)tC(3)*DDC(3),A(3*3)*A$(3,3) 

C 

c cohputfs the contribution of each point to thf observation equations 
c 

DXsXO-XC 

DYsTO-YC 

D7=?0“7C 

DR (1,1 )s-A(l, 1 )*DX"A(1*?)*0Y P A(1,3)#D7 
Dp(?,l)® - A(2,t)*DK»A(2*2)*nY"Af2,3)*D7 
DR(3,l)=-A(3,l)*DX»A(3,2)*nY-A(3,3)*0Z 

DP ( 1 , 2 )s SC AL * A ( J , 1 ) 

DP ( 2# 2 ) s SCAL*A(2,i) 

OBf 3,2 )b SCAL*A ( 3, 1 ) 

DP( l » 3 )b SC AL * A ( 1 , 2 ) 

DB(?,3)s SCAL*A(2,2) 
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DB( 3,3). SCAL*A(3,2) 

DB(l*4)e SCAL*A(1,3) 

DB(2»4). SCAL*A{2,3) 

DB( 3 » 4 )» SCAL*A< 3 , 3 ) 

DBf 1 » 5 )» SCAL*<Af 1 » 3 )*DV-Af 1 # 2 )*DZ> 

DB{?, 5 >* SCAL*<A( 2 , 3 )*DY-A( 2 , 2 )* 0 Z) 

0 B( 3 , 5 )* SCAL*( A< 3 , 3 )*DY-A( 3 , 2 )*DZ> 

OPC 1>63«-SCa^-*CaSC !#J ) *P V + aSC 1 , 2 )*OY + *S C l , 3 >*DZ > 
DBC2,6).-SCAL*<AS(2,1 )*0**AS(2»2)*DY+A$(2,3)*DZ> 

DB( 3#6 5*-ScAL*( AS (3*1 ) *[>* ♦ AS< 3, 2 >*DY ♦ AS ( 3, 3 > *DZ > 

DB( l, 75 *-SCAt*( A( 2 ,l)*Dy*A< 2 , 25 * 0 Y+A( 2 , 3 )*DZ) 

0P(2,n* SCAt*(A(l#l)*Oy 4 Af l# 2 )* 0 Y 4 A(i* 3 )*DZ) 

0B(3»7)« 0.0 
DO 10 K*l,3 
DO 10 1.1,3 

10 0DB(K,l )*»$C AL* At K ,L ) 

DDC(l) e XO-XOO 
DPC ( 2 ) s YO-YUO 
DPC ( 3 ) * 70-Z00 

C(l) . SCAL*(A(1,1)*DX+A(1,2)*0Y*A(1,3)*DZ)-X 
C f 2 ) * SCAL*<A(2,n*DX + A<2»2>*pY + AC2,3)*0Z)-Y 
C C 3 ) * SCAL*(A(3,1 )*DX+At 3,2 )*DY*A( 3, 3 )*DZ)-Z 
return 
END 

SUBROUTINE OFTCAL (Np, j,DP,ODB,DC,ODr,C,Nj,DN»0DW,0N,DDN,0K#DDK,BN, 
1 DDNINV.DDNIN j) 

C 

* 


T H I$ SUBROUTINE computes The cONTrIpUTTON OF each POINT To THF NO- 
RMAL equations of the form 

DN * DD ♦ BN * DOD * DK 
BNT * DO ♦ DON * DOD s DDK 

m atrI c fs cO m f>ute d by Thu subroutine * dn,oo n »bn,ok,dck 

INPUTS 


DB*DDB»DDC,C * EiFhENTS of ThE OBSERVATION EQUATIONS 

COMPUTED in subroutine part 
dc - element of the observation Equations compu- 
tet IN SUBROUTINE 3 DTRAN 

WJ,DW,DDh - wtlCHT MATRICES OF THE MODEL COORDS, THE TRA- 
NSFORMATION parameters and The groupo coords 

J - NUMBER OF THE POINT 


OUTPUTS 

pn,ddn,dk»ddk - eiFments Of The normal Equations for the 

INDIVIDUAL POINT J 

bn - element of the normal Equations stored in a 

NfP*3*3 MATRIv for FUTURE USE 
c DDNINV - INVERSE OF DON STORED IN A NCP*3*3 MATRIX 

c ddninj - inverse of ddn oe the individual point j 

C 

* 

c 

DIMENSION PK(7),ON(7,7),PDC(3),DPN(3,3),DB(3,7),HJ(3,3},p,0(7,7), 

1 CC3),DC(7),D0B(3,3)»0DNIF Jf 3,3),DDWf 3,3) 

DIMENSION DBT<7,3),DDBT(3,3),0BtNJ(7,3),DDBTNJ(3»3),BrJ(7,3),DDK1( 
13),DDK2(3) 

DIMENSION DDNINV(100,3,3),HN(100,7,3),DOK(100,3) 

C 
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c computes the contribution or each point to thf normal foliations 
c 

c 

c 

C INITIALIZE ELEMENTS Of ON DON p DK p DDK 
C 

DO 10 Kslp7 
0K(K 5=0 ,0 

DO 10 1*1*7 
10 DN ( L ) c 0.0 
DO 20 Ksl,3 
DOK ( JpK )*0, 

OP 20 L s 1 p 3 
20 DON(KpL) s 0,0 
C 
C 

C CONSTRUCT MATRICES 0BT{7#3) AMU DOB T ( 3p 3 ) 

C 

C 

DO 30 K*l;3 
DO 30 I = 1p 3 

30 DOBT(KpL> e DDB(LpK) 

DO AO V~\ »1 
DO AO I s 1 p 3 
AO D P T < K p t ) s DB ( L p K ) 

C 

C COMPUTE ON ( 7 , 7 ) AND DPN(3p3) 

C 

CALL MXMULT(DBTpWJpDRTNJp7p3p3i 
CALL MVMULT(DBTNJ#0B,DN,7^7,3) 

C 

CALL MXMULT (DDBTpWJpDDBTMJp3p3p3 ) 

CALL MXMULT(DOBTWjf DDB»DPN» 3*3,3) 

DO 50 K = 1 p 3 
Op 50 L sip 3 

DON(Kpl ) * DDN(K»L) + OPK(K,l) 

50 DDNINJ(KpL) s DON C K , L ) 

C 

CALL IN'VERS C DDNINj,3pII 1 
C 

00 60 K = lp3 
00 60 LeW3 

6C DONINVC JpKpL) * DDNINJCKpLI 
C 

C COMPUTATION Of DK(7) ANO UfK(Jp3> 

C 

OP 70 K * 1 p 7 
DO 70 L s 1 p 3 

70 DK ( K ) s DK(K) ♦ OBTWJ(K,l )*C(L) 

OP 60 K s 1 p 3 
DDKl(K) » 0,0 
DO 60 L s 1 p 3 

60 DDKl(K) = DDKl(K) ♦ DD& T ►' Jf Kp L 3 * C ( L > 

DO 90 Ksl , 3 
DDK? ( K 3 * 0,0 
DO 90 L s 1 p 3 

90 DDK2(K)bDDK2(K)4 DOW ( K *1 3 *ODC { L 3 
DO 100 K B 1 p 3 

100 DDK(JpK) b DDKKK) • OOK2CK) 

COMPUTATION OF BN(Jp7p3) 
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CALL MXMULT ( DBTWJ#DDB#BNJ, 7 # 3 # 3 ) 

DO 110 ** 1#7 
DO 110 L*l >3 

110 BNtJ#K#L) » BNJ(K#L ) 

RETURN 

END 

SUBROUTINE SIGMA ( NP#XP#Yn#ZU#XOO#YOO#ZOO#XC#YC#ZC# StAL#0ELTA# 
1 DDELTA#A>AS#X#Y#7'0C#WJ#PDW#DW>SIGNFW ) 

c 

c .... 


c 

C THIS SUBROUTINE COMPUTES THE MEAN SrUApE ERROR Of UN II WEIGHT Of 

EACH ITERATION# CALLED Sir, NEW# BY SOLVING THE OBSERVATION EOllA- 

t ions eor the computation of the residuals. 


INPUTS 

NP(NCP) 
X0#Y0#Z0 
XOO# YOO# ZDO# 

XOO#YOO,zOO 


X# Y#Z 
SCAL 


XC#YC»ZC 
DMt#pHl#CApA 
A, AS 
DELTA t 7#l) 
DDFLTA(Ncp,3) 
CONTINUE 

DC 


Wj,Dw,DDw . 


THf NUMBER OF CDNTROI S 

IHF GROUND COPRDS Of THE CONTROLS 

approximations to thf ground coords of the 
controls as used in the first iteraiion 

MOFEl. COORDINATES OF THE POINTS 
IHF APPROXIMATE SCAlF OF THE MOOEl 
APPROXIMATE TRANSLATIONS OF THE MODEL 

approximate rotations Of the model 
T pF rotation matrix ano ITS DERIVATIVES 

TpF COMPUTED TRANSFORMATION pArAmETERS 
Tpp cOmpUT pp GROUND COORDS Op The CONTROLS 

IpF RIghT hANP rATrIx Op ThE Op$ERv*TInN 
EfSl'ATlONS ( FROM SUBROUTINE 3pT R AN > 

weight matrices Or the muoel coords. The tra. 
NSfORmATTON PARAMETERS And ThE GROUND COORDS 


OUTPUTS 


signew - the ne w mean sou*he error of unit weight of 

the SPECIFIC ITERATION 


DIMENSION DBC 3 # 7 )#DDB( 3 # 3 )# ct 35 #DDCF 3 ),A( 3 # 3 ) # A SC 3 # 3 ># Dwl 7 # 7 ) 
DIMENSION DELtAF 7 )# DCf 7 )#V( 3 ) » DDV ( 3 ) # DV ( 7 ) 

DIMENSION X00(lO0)#X0(lOO>,Y0TlO0)#7OrlO0 >#X(100 )»Y(100)#Z(l00)# 
1 YOO(lOO)»70D<lOOl#WJf 100#3»3)#DDWC100#3#3)#ODELTA(100#35 
C 
C 

VTWVsO.O 
NDF * 3*NP*7 
C 
C 

J*l 

999 CALL PARK A # AS# XO ( J ) » YO ( J ) # ZD ( J ) # YOO ( J )# YOD( J ) # ZOD ( J ) # XC # YC # ZC # 
1SCAL#0P#DDB#C#DDC#X{ J)#YfJ)#Z( J>) 

DO 100 K*1 # 3 
Sl* 0 ,0 
S?=0,0 

DO 110 L=l#7 

110 S1»SW08(K#L)*DELTA(L) 

DO 120 L*1 # 3 

120 S2*S2 + D0B(K#L)*D0ELTA( J#L ) 

V ( K ) sC C K )*5 1 *S2 
100 DDV(K)*DDC(K)+DDELTA(U»K) 



o r> o o 
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SleO.O 
S2e0 ,0 

pn 130 k=i# 3 
DP 130 L*l#3 

SlsS1+V(K)*WjCj,K,L)*V(l 5 

130 S?=S?+POV(K)*DDW( J# K»L)*rOV(L) 

VTWV*VTWV*S1*S2 

lE(J-NP) 1 31 # 1 3?/ 1 3? 

131 J=J+1 

GO TO 999 
13? DP 190 K*1 * 7 
190 OVCKlsPCCKl+DElTACK ) 

DO 191 K ® 1 * 7 
DC Ml L e 1> 7 

14! VTNV*VTHV + PVfK)*PW(IOL)*rV(L) 

SIGNFW s VTWV/NDF 

RETURN 

END 

surroutine ROTATE (OHt'FHl »CApA, A, AS ) 

c 

C COMPUTES ThE rotation M A TklY FOR gJvEn RnTATInN ANqUE-S AND T pE 
C PApTIAt DERIVATIVES OF ThF FOTATION MATRly ELEMENTS FOP FHl 
C 

DIMENSION AC3#3 )»AS(3#3) 

SIND e SlN(OME) 

cost = cus(ome) 

S I Np s SIN(PKI) 

en$p & cuscphi > 

SINC = SlNfCAPA) 

COSC e COS(CAPA) 

AC 1 s> 1 ) = CDSP*COSC 

A ( 1 p 2 ) s CD50*SINC4-SIN0 *SINp*CUSC 

A(l>3>* SlNO*SlNC-COSO*SlNP*cnsC 

A ( 2 # 1 ) e *COSP*SINC 

A ( 2 p 2 )* C 0 $ 0 *C 0 SC-SIN 0 *SINp*slNr 

A ( 2 > 3 ) c SIN0*C0SC+CDSU*SINR*SINC 

A{3/.l)e S I NP 

A C 3* 2 ) = *-SIND*COSP 

A( 3 » 3 )s COSO*CO$P 

AS(1»!) ■ -SINP*COSC 

A S ( 1 # 2 ) a SINO*CPSP*COSC 

A 5 ( 1 » 3) e -CDS0*CDSP*CUSC 

AS( 2 * 1 ) " SINP*SINC 

AS(2»?) * -SIN0*C0SP*S1FC 

AS( 2 # 3 l * COSO*CC$P*$INC 

AS { 3 # 1 ) * COSP 

A S( 3 # 2 ) s SIN 0 *SINP 

AS( 3# 3) * "COSOiSINp 

RETURN 

END 

SUBROUTINE OUTpUF f I TEp» A S I GN , DElT A # SC AL» XC » YC* ZC # Opt * PHI » CApA, DE 
1 VPAR,SIGOEL # XO, Y0 # 7O*D0El TA,NUM pNCP#SIGMj) 

C 



c 

c this subroutine prints cut the resuits of the transformation 

c ApTFR ThE SOLUTION hAS CPNvErgEo.. IT ALSO PRINTS OUT ThE ADJUSTED 

values of the ground coordinated of the control points and their 

DEVIATIONS, 

Mllff Him 

c 
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C 

C 


DIMENSION 

DIMENSION 

dimension 


DELTA(7)#DEVHAP(7,7),SIG0FL(7*7> 

XO(100)#YO(IQO)#20(100>#DDELTA( 100#3>#DEVIA(100#3*3) 
SIGMJ(100#3# 3)#NUM( 100) 


NRITE{6,480) 1TEP 
480 FORMAT("l"#" ****** 

t DM NUMPER"#I3#/#12X 
1 ) 
WPITE(6»«90) ASIGN 
490 FORMAT (////#20X#" * 
1,** *",/#?0X#" *"#3X 


*********************** flNAL OllTPliT ********** 

-#//#i 2 x," the solution converges a r tep ite*ati 


#"%1gMa*ZeRD«"#F!1%'3X," *"#/#2 OX #"**"# 29X#"* 


^SCAL * SCAL ♦ DFLTA(l) 


AXC i XC ♦ DELTA ( 2 ) 

A VC * YC ♦ DELTAC3) 

A7C ■ 7C ♦ DELTAT4) 

AOME * OME * DELT A ( 5 ) 

ApHI * pH I ♦ DELTAC6) 

ACAPA e CAPA ♦ DFLTA(7) 

CALL RTODmS(GmE , I DEOmE# M l NOME# SECOmE ) 

CALL RTODmS(PHI , I DEPH I #M 1 NpHI # SECPH I ) 

CALL RT0DM5(CApA, I DC A P A , Ml C ApA * SF C APA ) 

CALL RTDDMS(DELTA(5)*IUrOME*MUnME #«DOmE ) 

call rtodmS(Delta(M, ioopht,mdphi #fdphi ) 

CALL RTODM$(DELTA(7)*lODrpA,MDCApA*SDCApA) 
CALL RTOUMS(ADmE , I D A 0 ME# M ADmE #8A0mE ) 
call Rtodmscaphi , idaphj,m#phi #saphi ) 

CALL RT0DMS/AcApA,lDACFA*MACApA,5ACApA ? 

CALL RTODMS(DE v pARc5,5)*lDOEvO*MlDEvO,SEDE v O) 
CaLL RT0l>MS(DEVPARt6*6 J# iO^EVP.MtDeVP^eDeVP) 
CALL RT0DMS ( DEVPAR(7*7)#I00EVC,MIDEVC,SEDEVC) 


500 

510 
520 
530 
5 AO 


1 


RAN SFOr HATI d N PARAMETERS"* 



DEVIATION",/) 


WPlTE(6#500) 

FORMat</'10X>" t 

"*nin VALUE"' 3 *'" C0 RR e C ti °N"»2Xj»" FINAL VAluE"*3x," 

WRITE (6, 510 )SCAL* DELTA (! )» ASCAL'OEVPARf 1# 1 ) 

FPRMATT5X," SCALF F AC TOR" # A X* 2 ( 2X , F) 1 ,5 # 2X , E 1 1 . A ) ) 

WRITE (6,520)XC ,DELT A ( 2 ) # AXC # DE VPAP (2 ,2 ) 

F0PMATT5X'" TRANSLATION IN X">2C?X»Fll»3'2X»Ell 
HR I TE(S»530)YC »DELTAC3)' AYC#DEVPAP(3#3) 

FORMAT (5X# " TRANSLATION IN Y",?{?X#FU«3#2X#E11,4)) 

WRITE<6,5A0)ZC ,0ELTAU>'AzC*0EVPAPU,4) 

F0RMATT5X*" TRANSLATION IN Z"#?(2X#FlI,3#2X#E11.4)) 

WRITE (6,550 > 1DE0ME,MIN0ME,SEC0ME# lrDOME,MOOME 
1 MAOmE ,SAOME*IDDEvO,MIDFvO»SEOEvO 
>50 E0RMATT5X*" ANGLE OMEGA", 5X, 4 f 1 X# 1 3# I 3' F6 , 2 ) ) 

WR I TE ( f , 56® ) IDEPHI ,MINPHI,SECPHI ,IOOpHl,MOpHI 
1 MAPHT ,SAPHI* I DDE VP# MI UE VP# SE DEV P 
>60 FORMAT ( 5X* " ANGLE PHI"' 7 X, 4 ( 1 X # I 3 » I 3# F6 « 2 ) ) 

WRITE (6,570) I DC A PA ,M1 CA PA ,SE C A pA , I PDCPA , MOCA PA ,SDCApA , I DA CPA , 

imacapa,sacapa#iddevc#midevc,sede VC 
>70 FORMAT (5 X' " ANGLE CAPA"# 6 x# 4 {1 x# I 3# I 3# F6 , 2 ) ) 

WRITE(A#580> 

>80 FORM AT < ///#10X#" THE VARIANCE-COVARIANCE MARTIX OF tml transformat 


.4)) 


.4 )) 


.A)) 


# SD0ME 

# idaome 

#SDpHI 

,I0APHI 

#SDCAPA 

,IDACPA 


1 IDN* # // ) 

DO 590 K«l*7 


590 WR 1 TF ( 6# 600 ) ( SIGOEL(K#l >'L*1#7) 

600 FORMATT" "'7EU.4 > 

WPITE(6,610) 

610 fORMAT(" 1"#//#18X#" LIST OF ADJUSTED GROUND COOROlNATES"'/# 18x* 


# m 
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.......................... NUMBER"* 6X*" X"*I1X*" V**, 

!1!X>" 7"*6X*" SIGCX)"*4X*" SIG(Y)"*AX," SIG(Z)"*//} 

DO 670 J=1 * NCP 

XO(JT » XO(J) + ODElTACj#l> 

YD(J) = YO(j) * DDF LT A ( J* 2 ) 

?D(J> = 70( J ) ♦ DDFLTACJ*3) 

DP 675 K * 1 # 3 

675 DEVIA(J,K*K> = SORT(SIGNJ( J#K*X)) 

6?0 WP I TE ( 6 * 6 30 ) NUM { J ) , XO ( J ) * Y0( J ) * 70 ( J > # DlV I A ( J * 1 * 1 >»DEVIA(J*2»?)*DF 
1 V r A c J*3*3) 

6 30 F0*HATH5*3(?X#Ml*3)»3(?X*EU*O) 

650 CONTINUE 
RETURN 
END 


C 

C 

C 

r. 

c 

c 

c 

c 

c 


SUPROUTINE OUTPUT ( I TE P , AS I GN , DFLT A , SC AL , XC , YC ,ZC , OmE , PH I ,C ApA , 
1DFVPAR ) 


this subroutine prints put the rfsuiTs of each iteration *xcEpT 
the last and final one* 


DIMENSION 0ELTA(7)f DEVKAPC7,7) 
WRITE (6* TO) ITER 

10 FORMAT f " 1 "*" ***************** 


ITERATION NUMBER"* J4,2x* 


1 " 

WPITF<6*?0) ASIGN 
70 EqRMATT ////*20X*" 


>*/>?0X*" 


*"*29x 


1*" *"*/*?0X*" *"*3X*" S I PM A ZEROe",F! 1 ,4#3X," *"*/,20X*" *"*?9X*" 


ASCAL s SCAL + DFLT A ( I ) 
AYC » XC ♦ DELTA ( 2 > 


AYC o YC ♦ DELTA f 3 > 


A?C u ZC * DELTAS) 

APMF * OME + DELTAT55 
APHI b PHI * DELTA(6) 

ACAPA c CAPA + DELTA<7> 

CALL RT0PmS( 0ME * IDEOME* NlNOME* SFCO mE ) 

call rtodms(Phi *idephi*njnphi,secphi) 

CALL RTOUMS(CApA,IDCApA,HlCAPA,SFCApA ) 

CALL RT00MS{DELTA(5)*I0rEMF # MDPMF *SPDmE ) 

Call RT0DMS(DELTA(6)*|0rPHl,M0PHT »SDPHI ) 

CALL RT00MS(DELTA(7),IDDCPA # mDCAPA,S0CAHA) 

CALL RTOOHS(AOmE * I D A 0 MF * M A OME *$AOmE ) 

CALL R700MS(APHI * 10 APHI *MAPM1 *SAPH] ) 

CALL RT0DMS(ACAPA,I0ACPA*MACAPA,SACAPA1 
CALL RTOPMSCDEVPARC5»5)*IDDEVO*MIDEVO,SEDEVO) 

CALL RT0UMS(DEVPAR(6#6),!PnEVP*MlDFVP*SEDEVP) 

CALL PT0PMS(DEVPAR(7»7 3* IDUE VC*MIDEVC , SEDEVC ) 

WPITE(6,30) 

30 F0RMAT(/#10x*" TRANSFORMATION PARAmF TEp S"» 

1/>10X," *//*25X, 

"" OLD VALUE", 3X," C ORRE C T I ON " , 3X * *' NEW VALUE"*3x»" DEVIATION",/) 
WRITEC6*40) SCAL*DELTA(l3*ASCAL*DEVPARf 1*1) 

AO FfiRMATfSX*" SCaLE F AC TOP"* AX* 2f 2X* E 1 1 , 5* 2X» E 1 1 . 4 ) > 

WPITE(6,50) XC *DELTA(2)#AXC #DEVPARf2*2) 

50 F0RMAT(5X*" TRANSLATION IN X "* 7 ( 2X , Fl 1 . 3* 2X * E V 1 , 4 ) ) 

WP I TE ( 6 * 60 ) YC *DELTA(3)* AYC ,PEVPAR(3*3) 



r> n n n n n on n r> 


c 

c 

c 

c 

c 
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60 F0RMATC5X#* TRANSLATION IN Y w # 2 ( 2X# F 1 1 , 3/ 2X# El 1 . A ) ) 

WPITF<6#70) ZC #DFLTA( 4 )# AzC # 0E VP ARf 4# 4 3 
70 FORMAT ( 5X# " TRANSLATION IN Z"# 2 C ?X # F 1 1 . 3# 2X # E tl , A ) ) 

WR I TF (6 # 80) IDEOME#MlNrME»SECOMF#IDDOME#MPOME #SOOMF #IDAOME# 

IMAOmF # $aome# iddfvo#midfvo»sfdfvo 

80 F0RMATC5X#" ANGLE 0MFGA*#5X# 4< IX# I )»I 3# F6.2) ) 

WRITF(6, 90) IDEPHI,MlNpHI,SECPHI#IPOPHl#MDPHI ,SDPHI ,II)APHI # 
1MAPHT ,SAPHI#ID0EVP»MIDFVP>SEDFVP 
90 FORMAT ( 5X# ** ANGLE PHI"# 7* # 4(1 X # I 3# I 3# F6 * 2 ) ) 

WRITE (f ,100) IDCApA,MlCAPA,SECApA,lrDCPA,MOeApA,SOCApA#IDACPA, 

1 MAC A PA# SAC A PA# I DDE VC# MlDFVC# SEOEVC 
100 F0RMAT(5X#» ANGLF CAPA"# 6X,4(1X#I3#I3#F6,2)) 

WPITE(6#1 10) 

no FORMAT f//#5x* " ROTATIONS IN radians*#//) 

WRITF(6,120) owe #0ElTA(5 AomE»DEVPAR(5#5) 

120 FPRMAT(5X#" ANGLE OMF G A* # 5X# 2 ( ?x# F II • 8# 2X# F 1 1 . A ) ) 

WRITF(6#130) PHI#DELTACf )#APHI#0FVPAR(6#6) 

130 FpRMAT(5X#* ANGLF PHI " # 5* # 2< 2X #F 11 . B# 2X # F 11 « 4 ) ) 

WRITF(6#1A0) CAPA#DFLTA(7)#ACAPA#DEVPAP(7#7) 

HO F0RMATC5X#* ANGLF CAPA " # 5x # 2{ 2X # F 1 1 . 8 # 2X# F 1 1 , 4 ) ) 

RFTURN 

END 

SURROUT I NE OUTPUX(JJ#NP# K CP#MA)(ITR#JPRINT#IPRINT#SIGO#OSIG) 

THIS SUPROUTINE PRINTS TrF INPUT OF THE MAlN/THREtD 
TEMPORARY WRITE-UP 
DIMENSION 0SIG(7#7) 

WPITF(6#10) NCP#MAXITR 

10 F ORM A Tf «l"#10x#* g^NErAi 0 & TA«,//10y,» NUMpER Or pOlNlS USED IN Tr 
IF SOLUTION*, 10x,l5#/,10y # * MAXIMUM NUMBER OF ITERATIONS AlLOWFD",! 
.1 1 X # 1 5 ) 

no forIaT(////2 3 x#* the v a p ian cE-covartancE matrix eor the hansforma 

1TI0N pArAmFTErS"#//# 10X#" (ORDER OF VARIABLES s SC ALE # OX# P V# D7# OmF 
l#PHI#CApA>",///> 

00 170 K* 1 # 7 

170 WRITF(6#1B0) (DSIG(K#L)#t *t#7) 

180 FORMAT (?X#7E10»3) 

ASIGO « SORT ( SI GO ) 

WRITE(8,110) ASIGO 

no FORMAT (////#20X#* ********************************"#/#20X>" 

),* *"#/#20X»* *"#4X," SIGMA ZER0*>"#F9,5#4X»" *"#/#20X#" 


* "# 29x 
*"#?9x#" * 


C 

C 


RFTURN 

END 


SUBROUTINE pANDU <1 X # I Y # YFL ) 



pANO 
RAND 
, R AND 

540 

10 

20 




rand 

30 

SUBROUTINE RANDU 



RAND 

40 




PANO 

50 

PURPOSE 



PAND 

60 

COMPUTES UNIFORMLY DISTRIBUTED 

RANDOM REAL 

numbers between 

RAND 

70 

0 AND 1.0 ANO RANDOM INTFGERS 

bftween zero 

AND 

PAND 

80 

?**31, EACH ENTRY USES AS INPUT AN INTEGER 

RANOLM NUMBER 

RAND 

90 

and produces a new intfgfr and 

RFAL RANDOM 

NUMpER, 

PANO 

100 




PAND 

110 

USAGE 



PANO 

120 
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call randu<ix,iy*yfl> 


description or PARAMFTERS PANP 150 

Ty - FOR THE FIRST ENTRY THIS MUST CONTAIN ANY EDO INTEGER RAND 160 
NUMBER with NINE OR LESS DIGITS, after the FIRST EnTr Y >PAND 170 
IV SHOULD bE THE PREVIOUS v*LUE OF lY COMPUTED RY THIS PAND IPO 
SUBROUTINE, PAND 190 

lY - A RESULTANT INTFGER RANOCM NUMBER RE 00 1 RE D FOR THF NfXTBAND ?00 
ENTRY TO THIS SUBROUTINE. THE RANGE OF THIS NUMpF R IS pAND ?l0 


CONTINUE 


BETWEFN 7ERU AND 2**31 

vfL- ThE RFSULTANT UNlpOPMt Y DISTRIBUTED# FLOATING POINT, 
'RANDOM NUMBEP IN THE RANGE 0 TO 1,0 


PAND 220 
PAND ? 30 
PAND ?40 
RAND 250 
PAND 260 
RAND 270 


REMARKS PAND 260 

THIS SUBROUTINE Is SPECIFIC to SySTEm/360 AND WILL PRODUCE rand 270 
?**29 TERmS BEFORE REPEATING. ThE RFFERENcE RFLOW pi ScilSSFSRANo ?flO 
SEEDS (65539 HERE), RUN PRDBLFMS» AND PROBLEMS CONCERNING PAND ?90 
PANOOm DIgHS llsUr This gENFPATTON SchEmE, hAcLAREN And PAND 300 
MARSAGLIA, JACM 1 7 * p, 83*P9» DISCUSS CONGRUFNTJAL PAND 3lO 

GENERATION METHODS and tests. ThE USE Of Two GENERATORS Of pAND 3?0 
ThE PANDU T y pF, OFF FILLING a TAbLE AND ONE P lc K lNG fROm THrPAND 330 
TAPLE. Is Of rENErTT in some CASES. 655A9 HAS bEEN PAND 3«0 

suggested as a $efd whUh has better statistical properties pand 35o 

for HIGH order BITS op ThF GENERATED DEVIATE. pand 360 

seeds should re chdsen in accordance with the discussion pand 370 

GIVEN IN The REFERENCE BELOW, ALSO* IT SHOULD fct NOTED ThATRAND 380 


If FLOATING POINT PANDOM NUMBERS Arp DESIRED, as are 
A vAlLABLF from RAFDU, ThF RANDOM CHARACTERISTICS OF THE 
FLOATING POINT DL v T A TE 5 APF MODIFIED AND IN pAcT THESE 
CONTINUE 

DEVIATES HAVE HIGH PROBABILITY OF HAVING A TRAILING LOW 
ORDER ZERO BlT IN THEIR FRACTIONAL part, 

SUBFOUTINFS AND FUNCTION SUBPROGRAMS REoulREp 
NONE 

METHOD 

POWER RE S I OUE M£lF J OD DISCUSSET IN IBM MANUAL C20“8011> 
RANDOM NlJMpEP GENERATION AND TESTING 


IYbJ**65539 

IE(IY)5,6,6 

5 IYbIY*?147483647*1 

6 YEL=!Y 

YEl*YEi * » 465661 3F IB 9 

RETURN 

END 

SUBROUTINE GAUSS(IX,S#Am, V) 


RAND 390 
BAND 400 
PAND A 1 0 

PAND A 2Q 
PAND 430 
RAND 440 
PAND 450 
PAND 460 
BAND 470 
RAND 480 
RAND 490 
BAND 500 

PAND 510 
i RAND 520 

RAND 530 
PAND 550 
PAND 560 
PAND 570 
RAND 580 
PAND 590 


GAUS 390 
GAUS 10 


, GAUS 


SUBROUTINE GAUSS J’J us A q 

PURPOSE gjug 

COMPUTES A NORMALLY DISTRIBUTED RANDOM NUMBER WITH A GIVEN GAUS 70 

MEAN AND STANDARD DEVIATION GA US 80 

nSACr GAUS 90 

USA ^. , . GAUS 100 

CALL G A U S S ( I X , S . A M « V ) GAUS 110 



o o fry 
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C 

C 

c 

c 


c 

c 
c 
c 
c 
c 
c 
c 
c 

c 
c 
c 
c 
c 
c 
c 
c 

c 
c 
c 

AeO » 0 

Drt 50 I = l» 12 

CALL RANDU(Ix'Iy>Y) 

ix=Iy 

50 A* A+ Y 

V*(A-6.0)*S+AM 

RFTUrn 

END 

SUBROUTINE INvER5<**M>n j 
C 

C II c 1 OIVIDION BY ZERO 

c this routine computes thl inverse of an m x m matrix* 
c And destroys the original matrix. 

c A 7 FPO DIAGONAL ElEmENT Wilt CAUSE A DIVISION by ZERO 
c WHICH WlIL STOP the PROGRAM, 

C 

DIMENSION A(M»M) 

II e 0 

40 Do 1 L » 1 * M 
C ROW L IS PIVOTAL ROW 

IF(A(L*1) .LT. 1.0E-15) GO TO 8 
FMNs 1 . 0/A( L* I ) 

DlV = A(L,1) 

DO 4 J*2»M 

4 A C L > J-l>* A(L*J)/DIV 
A(L#M)b FMn 

c PIVOTAL row HAS pEEN MODlFlEr * NOW OTHER ROWS TO BE MODIFIED 
DO 7 J*l' M 
I F ( J-L ) 5# 7*5 

5 FMULT c A(J,1) 

DO 6 K*2> M 

6 A(J# K-l) «CACJ#K)-(CFMl'l T)*CA(L*K-im> 

A ( J» M ) b -FMN * FMULT 

7 CONTINUE 


DESCRIPTION OF PARAMETERS 

!X -IX MUST CONTAIN An ODD INTEgFR NUMBER WITH NINE OR 
LESS DIGITS ON THE FIRST FNTRY TO GAUSS. ThEREAFTFR 
IT Hi LI CONTAIN A UNIFORM) Y DISTRIBUTED INTEGER RANDOM 
number generated by thf subroutine FOR USF UN THE nfxt 
entry to the subroutine, 

5 -the dEsireo standard deviation of the normal 

CONTINUE 

distribution. 

a m -the desired mfan or the normal distribution 
v -the value of the computed normal random variable 

remarks 

this subroutine uses ranou which Is machine specific 

SUBROUTINES and FUNCTION SUBPROGRAMS REQUIRED 
RANOU 

mfthod 

USES 12 UNIFORM RANon M NUMBERS TO COmpUTF NORmAL RANDOM 
numbers by central limit theorem, thf result is then 
adjusted to match thf given mean and standard deviation, 

ThE uniform RANDOM NUMBERS COMPUTED wITHIN The SUBROUTINE 

arf found by the power rfsiduf mfthdd. 


GAUS 120 
GAUS 1 30 
GAUS 1 40 
GAUS 150 
GAUS 160 
GAUS J 70 
GAUS 180 
GAUS 190 

GAUS 200 
GAUS 210 
C,AUS 220 
GAUS 230 
GAUS 240 
GAUS 250 
GAUS 260 
GAUS 270 
GAUS 280 
GAUS 290 
GAUS 300 
C,AUS 310 
GAUS 320 
GAUS 330 
GAUS 340 
GAUS 350 
GAUS 360 
GAUS 370 
GAUS 380 
GAUS 400 
GAUS 410 
GAUS 420 
GAUS 430 
GAUS 44 O 
GAUS 450 


INVS 10 


I NVS 20 
INVS 30 
JNVS 40 

INVS 50 
JNVS 60 
INVS 70 
INVS 80 
INVS 90 
INVS 100 

INVS 110 
INVS 120 
INVS 130 
INVS 140 
!NV$ 150 
INVS 160 
INVS 170 
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RETURN 
ft 11 = 1 

return 

END 

SUPRDUT I NE RTOOMS ( A N G ( F , 1 DEG# MIN# SEC ) 

C 

C PURPOSE 

c transforms angles from radians td degrfes#minutes* seconds 
c 

c description of variables 
c angle - THE INPUT Angle td transformed 
c IpEfi - THE DEGREES f INTFgER ) 

C min - THE MINUTES ( INTEGER )) 

c sec - the seconds ( reaj > 

c 

ANG = angle* { 57,?9577VM3J ) 

IDEG * I NT( ANG) 

DEG = IDEG 

EMIN = ( ANG - DEG )* 60,0 
M I Ne I NT < EMIN) 

TMlN = Ml N 

SFC= ( FMIN - TM I N ) * 60.0 

rfturn 

END 

SUBROUTINE MXMULT ( A#B*C»K,L*M ) 

C 

C 

c 

C ThF SUPROUTINE MULTIPLIES t*o MATRICES 

C A PRF MULTI PLIER MATRIX WTTH DIMENSIONS K BY M 

C R POSTMULTIPLIER MATRIX WITH DIMENSIONS M BY L 

C C RESULT MATRIX WITH DIMENSIONS K RY L 

C 

C 

c 

DIMENSION ACK»M),B(M,l),C(K,L) 

C 

c 

DO 10 IrliK 
DO 10 Js 1 ; L 
C(I#J) =0,0 
DO 10 N*1#M 

10 C(Ij>J) = C(Ij>J) ♦ A{ I,N)*B( N# J) 

RETURN 

END 

SUBROUTINE PTCOMP ( Np #Nf P> SC AL * XC * VC > ?C » 0ME » PHI t C AP A# X# Y * Z* NUM* 
1 DEVPAR>SIGJ) 

C 

C , , 


this subroutine will compute the transformed coords for the 

UNkNOwN POINTS AFTEr THE ADJUSTMENT Of THE TRANSFORMATION 


DIMENSION NUM(100)*X<100)» Y{100)» Z{100)#A(3»3)j.ASC3*3) 

1 #DEVPAR(7#7)^ SIGj(100#3*3) 

Dimension ATC3»3)#XMC0OPMO0*3)#GRCn0Rf 100* 3>* ATX(100#3>,f R ANSL(3> 

C COMPUTE AND TRANSPOSE I HE ROTATION MATRIX 

C 

CALL ROTATE(OME#PHl»CAHA# A, AS) 


INVS 1 Bo 
INVS 190 


RTUMS 10 
RTDMS 20 
PTDMS 30 
RTDMS 40 
PTOMS 50 
PTDMS 60 
RTOMS 70 
RTOMS 80 
RTOMS 90 
RT0MS1 00 
PTDMS1 10 
RTDMS1 20 
PTDMSI 30 
RT0MS140 
PTDMSI 50 
PTDMSI 60 
PTDMSI 7 0 
RTDMS10O 
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DD 10 1*1*3 
00 10 J*l*3 
10 AT ( I * J) * A ( J* I ) 

C 

c FORMULATE model coordinate and TRANSLATION VECTORS 
c 

DO 20 JsNCP* NP 
XMC00R(J*1> * X( J) 

XMC00R(J*2) * YCJ) 

20 XMCPORf J* 3 ) * 7 C j ) 

TRANSL(l) * XC 

TRAN$Lf2) * YC 

TRARSLY3) ■ ZC 

ROTATE MODEL 

DO 30 JeNCP* NP 
DO 30 K*l *3 
ATX(J*K) « 0, 

00 30 L *1 * 3 

30 ATX(J#K) ■ AT X ( J* K ) ♦ AT ( K * L ) * XMC UOP< J * L ) 

compute ground COORDINATES of thf unknown POINTS 

DO 40 J * NCP* NP 
DO 40 I *1*3 

40 GRC OOR f J* I ) c (l/SCAL)*ATXf J* I ) + TPANRL(I) 

print t h e ground coordinates 0r the unknown points 

VLb DEVPAR(1*1) *DEvPARc 1*1) 
vXTs DE yRA R (2*2) *DFvpApf2>2) 

VYT*0EVPAR(3,3) *DEVPAR(3*3) 

V?T* DEvPARf4>4) *DFVPAPf4*4) 

V0bDfVPaR<5>5) *DeVP a R{?*5!i 
VP*DFVPAR c 6 * 6 )* 0EVPAR(6*6) 

VC° OEvpAR t 7* T ) *DE V PAH(T*T) 

WRITE (6*50) 

50 fORMAT^I "*////*10x*" LIST OF COMPUTE D GROUND COORDINATES OF L'NkNO 

l WN POINTS"*/* 10X*" 

NUMBER" * 1 OX #" X"*1*X," Y"*15X*" Z"* 

1 1 0 X» "SIGMA X"* 1 0 X* "SIGMA Y»* \ OX* "STGMA Z"* /) 

00 80 J * NCP*NP 
XMs XMC(jOR<j*l) 

YM*tXMCPOR< J*2) 

ZM»XMC00R( J*3) 

VXM«SIGJ( J, 1*1) 

V YMs SIGJ(J,2*2) 

VZM* SIGJ(J*3*3) 

VX*( (A (1 * 1 )*XM ♦ A(2*1)*YM ♦ A ( 3 * 1 )* Zm )**2 ) *VL /(SCAl**4) ♦ 

1 ( ( 1.0/SCAL)**?)*( ( ( A$f 1* 1 )*XM ♦ AS(?*1)*YM ♦AS(3»1)*2M)**2) 

1 *VP ♦ (<-A(2,1)*XM * A( 1 , 1 )*YM)**2)*VC +( A( 1 , J ) **? ) *VXM 

1 ♦ ( A(2# 1 )**2)*VYM ♦ ( A( 3* 1 )**?)*V2M ) ♦ VXT 

VXsSORT(VX) 

VY* (<A(1*2)*XM +AC?*2)*YM ♦ A ( 3 * 2 ) *Zm ) **2 ) * VL /(SCAL**«) 

1 + ( ( 1 « C/SC AL )**2) * 

1 (C(«A (1,3) *XM • A(2,3)*YM •A(3,3)*?M)**2)*V0 

1 ♦( ( AS( 1 »2)*XM -f AS(?*2)*YM + AS(3*2)*7M)**2)*VP 

1 ♦( ( At ?* 2)*XM • A( 1*2)*YM)**2)*VC 

1 ♦ (A(l*2)**2) * VXM ♦ (A(?*?)**2)*VYM ♦ ( A ( 3* 2 )**2 )* Y7M ) 

1 ♦ VYT 

VY« SOPTCVY) 
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V Z = (<A(1,3)*XM ♦A(2,3)*VM ♦ Af 3* 3 )*ZM)**2 )*Vt /(SCAL**A) 
t + <(1.G/SCAL>**2)* 

1 (<(A<1»2)*XN ♦ A(2»2)*YM + A ( 3 # ? )*7M ) **2) *V0 

1 +({A$(1,3)*XM +AS(2»3)*YM ♦ AS{ 3, 3 ) *ZM >**2 ) *VP 

t ♦({*<?, 3)*XM - A( 1#3)*YM)**2)*VC 

1 ♦ ( A( 1 ,3)**? )*VXM ♦ (A(?,3)**2)*VYM ♦ ( A ( 3 , 3 ) **? ) *VZH ) 

1 ♦ V7T 

V?sSQRT(V2) 

PO WRITF ( 6 # 90 ) NUM(J>»( G*C POP ( J , T > * I *1 $ 3 W X » V Y , VZ 
90 F0PMAT(5X#I5*6(6X*ntt3)) 

RETURN 

EK’D 
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6. SAMPLE RUNS 

6.1 TO PERFORM ABSOLUTE ORIENTATION 


6.1,1 INPUT DATA 


2 32 6 6 0 t 

.0001 

1.000E+20 1.000F+20 l.OOOE+^O 1.000F+20 1»OOOE*20 1.000E+20 



27 29 

15 

14 

27 

001806.0 

015096,0 

310,0 

29 

010759,0 

021624,0 

386,0 

15 

014080,0 

-009575,0 

904,0 

14 

008007,0 

-006116,0 

360,0 

17 

00«701.0 

005200.0 

301.0 

3o 

017232,0 

009647,0 

347,0 

27 

-045,406 

058,396 

-000.799 

29 

014,964 

100*865 

-000,061 

15 

033,930 

• 1 07 « 6 35 

003, 133 

14 

•006,402 

-083.653 

-000,215 

17 

•000,479 

-008,434 

-000,704 

30 

056.822 

020.231 

-000.362 

2 

049,594 

036,132 

003,494 

28 

-038,419 

106.586 

*000.420 

175 

-008,057 

-122.278 

000.275 

182 

019,160 

-117.189 

003,096 

190 

056,678 

-115.561 

004,227 

181 

019,586 

- OP 4 , 599 

001.958 

174 

-010.957 

•088,442 

•000.351 

166 

-042,525 

-086,624 

■000,229 

17? 

-010.956 

-021.763 

•000, 583 

173 

-013,266 

*057.254 

-000.367 

1 8 0 

020.133 

-050, 199 

000.056 

188 

057,172 

-055.845 

002,183 

169 

-009,653 

072,515 

-000. 170 

177 

020,824 

075,133 

-000.128 

184 

054,609 

076,321 

003,111 

185 

058,569 

047,865 

002,914 

178 

026,090 

052,083 

000,028 

170 

-010,685 

051.498 

-000,236 

163 

-041,656 

005,040 

-000,679 

171 

-005,245 

018,659 

-000, 126 

186 

054,842 

013.650 

-000,223 

187 

053.032 

•01«,I09 

000.507 

199 

028,979 

-023,733 

•001,182 

168 

-011.219 

108,819 

-000,041 

183 

058,296 

112,773 

002,585 

217 

023.897 

110,495 

000,339 

27. 

225E403 


, 1F+03 

29. 

225E+03 


.1F+03 

15. 

225E403 


• 1F*03 


1 . OOOE+20 i,oooe;*2o 


6,25E 4 Q2 

6.25E+02 

6.25E+02 
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H. 

225E+03 

, i r*o3 

6. 25E+02 


17. 

225E+03 

tlF+03 

6 * 251*02 


30. 

225E+03 

♦ 1F+03 

6.25E+02 

99 

27 

. 1 E“03 

. IF-03 

♦IE-03 


30 

• IE-03 

.ir-03 

♦IE-03 


29 

.IE-03 

•IF-03 

♦IF-03 


15 

.IE-03 

, 1F-03 

♦IE-03 


1 0 

.IE-03 

• 1 F*03 

♦ 1E*03 


17 

♦IE-03 

» IF-03 

♦IF-03 


28 

.IE-03 

.1F-03 

♦1F*03 


? 

.1F-03 

.IF-03 

• 1 E"03 


190 

.IE-03 

.IF-03 

♦IF-03 


162 

.IE-03 

♦IF-03 

♦IF-03 


1 66 

.IE-03 

.IF-03 

♦IF-03 


175 

• IF-03 

• IF-03 

♦IE-03 


17ft 

. 1F=03 

. 1F®03 

.1F-03 


161 

.11-03 

.IF-03 

.IE-03 


188 

.IE-03 

.IF-03 

♦IE-03 


180 

.1F-03 

. IF-03 

♦IE-03 


173 

.IE-03 

♦IF-03 

* 1 E — 0 3 


172 

.IE-03 

♦IF-03 

• 1 E"03 


199 

.1F-03 

♦IF-03 

♦IE-03 


167 

.IE-03 

♦IF-03 

.1E*03 


166 

.1F-03 

.IF-03 

♦ IE-03 


171 

. IE-03 

♦IF-03 

.IE-03 


163 

. IE-03 

♦IF-03 

♦IE-03 


170 

e SF-03 

.IF-03 

♦IE-03 


178 

.IE-03 

.IF-03 

.IF-03 


185 

•1F-03 

♦IF-03 

♦IE-03 


1 6ft 

.11*03 

. IF-03 

♦ 1 E*03 


177 

♦ 1 F"03 

♦IF-03 

♦ 1 F*03 


169 

• 1 E*03 

♦ 1 F"03 

♦ 1 E"03 


168 

• 1 F -0 3 

♦IF-03 

. IE-03 


183 

» 1 E*0 3 

♦lF-03 

♦IE-03 


217 

. lE*03 

♦ 1 F*03 

♦IE-03 

99 



6.1.2 OUTPUT . 

GFNERAL DATA 

NUMBER OF POINTS U&t T IN THE SOLUTION 6 

MAXIMUM NUMBER OF ITERATIONS ALLOWED 6 


thf variance -covariance matpix for the Transformation parameters 



(ORDER OF 

VARIABLES r 

SCALE, dx. 

py , nz , ome , 

PHI, CAPA) 


O.lOOF 

21 0,0 

0,0 

0,0 

0,0 

0.0 

0,0 

0,0 

O.lOOF 

21 0,0 

0,0 

0,0 

0,0 

0.0 

0,0 

0,0 

O.lOOE 21 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

O.lOOE 21 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

O.lOOE 21 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

O.lOOE 21 

0,0 

0,0 

0,0 

0.0 

0,0 

0,0 

0,0 

O.lOOF 21 



* sigma zero* q.ciooo * 
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********************* OUTPUT OF SUBROUTINE I NPUT2 ******************* 


T R A N 

S F 0 

R 

M A T I 0 

N P A R A M 

E T t R 

PARAMETER 




COMPUTED APPROXIMATIONS 





IN PAD 

IN D 

SCALE 




0.00667 

TRANSLATION IN 

X 


10066 

.594 

translation IN 

V 


5982 

.703 

TRANSLATION IN 

z 


434 

.501 

ROTATION 

omega 



0.00022 

0 

POTATION 

PHI 



0,00166 

0 

POTATION 

CAPA 



0,01614 

0 1 


M 0 

0 

EL CO 

n P 0 I N A T 

t s 

1. CONTROL 

POINTS 





NUMBER 

POINT 


X 

y 

7 

1 

27 


•45,4060 

58,3960 

* 0,7990 

? 

29 


14,9640 

100,8650 

"0.0610 

3 

15 


33.9300 

-107,6350 

3.1330 

4 

14 


=6,4020 

*63,6530 

■*0.2150 

5 

17 


-0,4790 

-8,4340 

*0,7040 

6 

30 


56,6220 

20.2310 

"0, 36?0 

2, UNKNOWN 

POINTS 





NUMBER 

POINT 


X 

Y 

Z 

7 

2 


49,5940 

36,1320 

3,4940 

6 

28 


-38,4190 

106,5860 

-0,4200 

9 

175 


-8, 0570 

-122.2780 

0.2750 

10 

162 


19,1600 

-117,1800 

3,0960 

11 

190 


56,6780 

-115,5610 

4,2270 

12 

161 


19.5860 

-84,5990 

1,9580 

13 

17 A 


-10,9570 

•88,4420 

*0,3510 

H 

166 


-42,5250 

-66,6240 

-0.2290 

15 

172 


=10,9560 

*21,7630 

-0.5830 

16 

173 


-13,2660 

*57,2540 

-0,3670 

17 

180 


20,1330 

*50. 1990 

0,0560 

10 

166 


57.1720 

-55.6450 

2,1830 

19 

169 


-9,6530 

72.5150 

-0.1700 

20 

177 


20.8240 

75,1330 

*0.1280 

21 

184 


54,6090 

76,3210 

3,1110 

22 

185 


5». *690 

47,6650 

2.9140 

23 

178 


26,0900 

52,0830 

0,0280 

24 

170 


-10,6850 

51.4980 

-0,2360 

25 

163 


•41 ,6560 

5,0400 

•0.6790 

26 

171 


•5,2450 

16.6590 

•0,1260 

27 

186 


54.P420 

13,6500 

-0,2230 


0 

5 

55 


46,253 
41 .500 
29,799 
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28 

167 

53.P320 

-19.1090 

0.5070 

29 

199 

26.9790 

*23.7330 

-1,1620 

30 

166 

-11.2190 

106,6190 

-0.0910 

31 

163 

56.2960 

112.7730 

2.5650 

3? 

217 

23.P97Q 

110.9950 

0.3390 
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GROUND COORDINATES 


CONTROL POINTS 


NUMBER 

POINT 

X 

Y 

7 

t 

27 

1 AO&.COOO 

1 5096.0000 

310,0000 

2 

29 

1075V, rooo 

21624,0000 

386,0000 

3 

15 

14080,0000 

*9575,0000 

904,0000 

4 

14 

8007,0000 

-6116,0000 

360,0000 

5 

17 

8701 * 0000 

5200, 0000 

301.0000 

6 

3 0 

17232.0000 

9647,0000 

347.0000 
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ThE ' 

VARIANCE 

•COVARIANCE 

MAIRI 

CES FOR 0 ROUND 

COORDINATE 

W 

f P 

Pfl I NT 

27 

0.2250Q0E 

03 

0,0 


0.0 




0,0 


0.100000F 

0 3 

0,0 




0,0 


0,0 


0 » 625000E 

03 

POINT 

29 

0.225000E 

03 

0.0 


n.o 




0,0 


0.100000P 

03 

0,0 




0,0 


0.0 


0.625000F 

03 

POINT 

15 

0.225000F 

03 

0,0 


0,0 




0,0 


o.iooooor 

03 

0,0 




0.0 


0,0 


0,625000? 

03 

POINT 

1 A 

0.225000? 

03 

0.0 


0,0 




0,0 


o.iooooor 

03 

0,0 




0.0 


0.0 


0,625000? 

03 

POINT 

17 

0.225000E 

03 

0,0 


0.0 




0,0 


o.iooooor 

03 

0,0 




0,0 


0.0 


0, 625000? 

03 

POINT 

30 

0,2250000 

03 

0,0 


0,0 




0,0 


o.iooooor 

03 

0,0 . 




0.0 


0.0 


0.625000F 

03 



the 

point 

PH I NT 
PHI NT 
PD I NT 
POINT 
POINT 
POINT 
POINT 
POINT 

point 

POINT 

POINT 

POINT 

POINT 
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V4RIANCE-C0VARIANCE MATRICES FOR MODEL COORDINATES 


27 

0.100000F-03 

0, 

0 

0,0 


0.0 

0. 

iooooor*o3 

0.0 


0.0 

0. 

0 

0. 100000E*03 

29 

0.1 OOOOOE*0 3 

0. 

0 

0.0 


0,0 

0. 

100000F-03 

0,0 


0,0 

0. 

0 

0.100000E-03 

15 

0,100000E*U3 

0, 

0 

0,0 


0,0 

0, 

100000F-03 

0.0 


0,0 

0, 

0 

0.100000E-03 

1 A 

o,iooooor*o3 

0, 

0 

0,0 


0.0 

0. 

100000F-03 

0,0 


0.0 

0, 

0 

0,1 OOOOOE-03 

17 

0.100000E-03 

0. 

0 

0,0 


0,0 

0. 

1000G0r-03 

0,0 


0,0 

0. 

0 

0 , lOOOOOE-o* 

30 

0, 100000E-03 

0. 

0 

0,0 


0,0 

0. 

I00000r-03 

0.0 


0.0 

0. 

0 

0, 100000F-03 

7 

0, 100000E-03 

0, 

0 

0,0 


0.0 

0. 

iooooor-03 

0.0 


0.0 

0. 

0 

0,1 000Q0E*03 

2fi 

0.1 00Q00E*Q3 

0, 

0 

0,0 


0,0 

0, 

1 OOOOOF'O 3 

0,0 


0.0 

0. 

0 

0. 100000E-03 

175 

0, IOOOOOE'03 

0. 

0 

0,0 


0,0 

0, 

iooooor-03 

0.0 


0.0 

0, 

0 

0.100000E-03 

1 82 

0.100000E-03 

0. 

0 

0,0 


0,0 

0, 

t00000r*03 

0.0 


0.0 

0, 

0 

0, 1 OOOOOE-03 

190 

0.100000E-03 

0. 

0 

0,0 


0,0 

0. 

i00000r»03 

0,0 


o 

m 

O 

0. 

0 

0.100000F*03 

181 

0.100000E-03 

0, 

0 

0,0 


0,0 

0, 

iooooor-o3 

0.0 


0.0 

0, 

0 

0.100000E-03 

17A 

0, 100000E-03 

0. 

0 

0,0 


0.0 

0, 

100000F-03 

0,0 


0,0 

0. 

0 

0.100000E-03 

166 

0, 100000E®03 

0. 

0 

0,0 


0.0 

0, 

o 

o 

o 

o 

o 

• 

o 

UJ 

0,0 


0,0 

0. 

0 

0.100000E-03 

172 

0.100000E-03 

0. 

0 

o 

* 

o 


POINT 
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0»0 

0,0 

o.looooor-o3 

0,0 

0,0 

o, ioooooe*o3 

POINT 

173 

0. 100000E-03 
0,0 
0,0 

0,0 

0, 1000007*03 
0,0 

0.0 

0,0 

0, 100000E-03 

POINT 

ieo 

0, 100000F-03 
0,0 
0,0 

0,0 

o , i ooooor*o 3 
0.0 

0.0 

0,0 

0. 100000E“03 

POINT 

tee 

0,100000E*03 

0,0 

0,0 

0,0 

0 , 1 OOOOOr “03 
0.0 

0.0 

0.0 

O.IOOOOOE-03 

POINT 

169 

0.100000E-03 

0.0 

0,0 

0,0 

0. 100000r-03 
0,0 

0,0 

0.0 

0.100000E-03 

POINT 

177 

0.100000E-03 

0,0 

0,0 

0,0 

o. iooooor-o3 
0.0 

0,0 

0,0 

0,i00000p"03 

Pn i nt 

180 

0.100000F-03 

0,0 

0,0 

0,0 

0.100000r-03 

0,0 

0,0 

0,0 

0, 100000F-03 

POINT 

189 

0.100000F-U3 

0,0 

0,0 

0,0 

0.100000r-03 

0*0 

0.0 

0.0 

0,1 00000E®03 

POINT 

178 

0. 100000F-03 

0.0 

0.0 

0,0 

o.iooooor-o3 

0.0 

0,0 

0.0 

0, 100000E»03 

POINT 

170 

0.100000F-03 

0,0 

0.0 

0,0 

0. 100000P-03 
0,0 

0,0 

0.0 

o.i ooooof*o3 

POINT 

163 

0.100000E-03 

0.0 

0,0 

oo o 

O o 

o 
o 
o 
o 
o 

9 

o 

'Jj 

0,0 

0,0 

0, 100000E-03 

POINT 

171 

0.100000t-03 

0,0 

0*0 

0,0 

o,iooooor-o3 

0,0 

0.0 

0,0 

0. 100000E*03 

POINT 

186 

0.100000F-03 

0,0 

0,0 

0,0 

0. 100000F-03 
0,0 

0.0 

0,0 

o.i oooooe*o3 

POINT 

187 

0,100000F*03 

0,0 

0.0 

0,0 

0.100000P-03 

0.0 

0,0 

0.0 

0, 100000F-03 

POINT 

199 

0.100000F-03 

0.0 

0.0 

0,0 

o.iooooor-oT 

0,0 

0.0 

0,0 

0.100000E-03 

POINT 

168 

0,100000F*03 

0,0 

0,0 

0, 100000F-03 

0,0 

0.0 





0,0 

72 

0,0 

0. 100000003 

POINT 

183 

0,100000003 

0.0 

0.0 

0.0 

0.100000003 

0.0 

0.0 

0.0 

0,100000003 

POINT 

217 

0.100000003 

0.0 

0.0 

0.0 

o,iooooooos 

0,0 

0,0 

0,0 

0 , 10000000 3 
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HEIGHTS USED IN SrlUTlON ********************** 


THr WFlGHT MATRIX FOR THf TPANSFORMAT fCN PARAMETERS 


*1 00F-23 

0.0 

0,0 

o.o 

0,0 

0.0 

0.0 

,P 

0, 100F-23 

0,0 

0,0 

0,0 

0.0 

0.0 

.0 

0.0 

0,100E“?3 

0,0 

0,0 

0.0 

0.0 

.0 

0.0 

0.0 

O.IOOE-23 

0,0 

0,0 

0,0 

.0 

0,0 

0,0 

0.0 

0.100E-23 

0.0 

0,0 

.0 

0,0 

0,0 

0,0 

0,0 

0, 100E-23 

0,0 

.0 

0.0 

0.0 

0,0 

0,0 

0,0 

0, IOOf-23 


The w fi ghT 

matrices pOR 

The ground 

COORDINATES Of TrE CONTROL POINTS 

POINT 

27 

0 t AA4AF-06 

0.0 

0,0 



0,0 

0»1POOE"05 

0.0 



0,0 

0*0 

0, 1600E"06 

POINT 

29 

0, 4444E-06 

0*0 

0.0 



0,0 

0» 1000E*05 

0,0 



0,0 

0,0 

0. 1600E*06 

POINT 

15 

0.4444E-06 

0.0 

0,0 



0,0 

0» 1 P00E*05 

0.0 



0,0 

0.0 

0, 1600E*06 

POINT 

1 A 

0, 4444E-06 

0,0 

0,0 



0,0 

0.1POOE-O5 

0.0 



0,0 

0.0 

0, 16O0E-O6 

POINT 

17 

0.4444E-06 

o.o 

0.0 



0,0 

0.1000E-05 

0,0 



0,0 

0.0 

0, 1600E-06 

POINT 

30 

0.4444F-06 

0.0 

0,0 



0,0 

0.1000E-05 

0,0 



0,0 

0.0 

0.1600E-06 


THF WEIGHT MATRICES FOR T MF COORDINATES tlF T HF MODDEL POINTS < CONTROL AND UNkN 
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POINT 

27 

O.IOOOE 

01 

0.0 


0.0 




0,0 


0 • 1 000£ 

01 

0.0 




0,0 


0.0 


O.IOOOE 

01 

POINT 

29 

O.IOOOE 

01 

o.r 


0.0 




0,0 


O.IOOOE 

01 

0.0 




0,0 


0*0 


O.IOOOE 

01 

POINT 

15 

0 . 1 OOOE 

01 

0.0 


0,0 




0,0 


O.IOOOE 

0! 

0,0 




0,0 


0.0 


O.IOOOE 

01 

POINT 

1A 

O.IOOOE 

01 

0.0 


0,0 




0,0 


0»1000E 

01 

0,0 




0,0 


0*0 


O.IOOOE 

01 

POINT 

17 

0 , 1 OOOE 

01 

0*0 


0,0 




0,0 


O.IOOOE 

01 

0,0 




0,0 


0*0 


O.IOOOE 

01 

POINT 

30 

0, 1000E 

01 

0.0 


0,0 




0.0 


O.IOOOE 

01 

0,0 




0,0 


o»c 


O.IOOOE 

01 

POINT 

7 

0,1 OOOE 

01 

0.0 


0,0 




0,0 


O.IOOOE 

01 

0.0 




0,0 


0*0 


O.IOOOE 

01 

POINT 

20 

0, 1000E 

01 

0.0 


0.0 




0,0 


O.IOOOE 

01 

0.0 




0.0 


o 

• 

o 


O.IOOOE 

01 

POINT 

175 

0 , 1 OOOE 

01 

0*0 


0,0 




0,0 


O.IOOOE 

01 

0.0 




0,0 


0.0 


O.IOOOE 

01 

POINT 

18 ? 

0 , 1 OOOE 

01 

0.0 


0,0 




0,0 


O.IOOOE 

01 

0.0 




0,0 


0.0 


o 

• 

o 

o 

o 

m 

01 

POINT 

190 

0 , 1 OOOE 

01 

0.0 


0,0 




0,0 


O.IOOOE 

01 

0,0 




0,0 


0.0 


O.IOOOE 

01 

POINT 

161 

0 , 1 OOOE 

01 

0.0 


0.0 




0,0 


O.IOOOE 

01 

0.0 




0,0 


o»o 


O.IOOOE 

01 

POINT 

1 7 A 

0,1 OOOE 

01 

0.0 


0.0 




0,0 


O.IOOOE 

01 

0.0 




0,0 


0.0 


O.IOOOE 

01 

POINT 

166 

0,1 OOOE 

01 

0,0 


o 

• 

o 




0,0 


O.IOOOE 

01 

0.0 




0,0 


0.0 


O.IOOOE 

01 

POINT 

172 

0,1 OOOE 

01 

0.0 


0,0 




0.0 


O.IOOOE 

01 

0,0 




0,0 


0.0 


O.IOOOE 

01 
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POINT 173 O.IOOOE 01 
0,0 
0,0 

POINT 160 0,1 OOOE 01 
0,0 
0,0 

POINT 186 O.IOOOE 01 
0,0 
0.0 

POINT 169 O.IOOOE 01 
0,0 
0,0 

POINT 177 O.IOOOE 01 
0,0 
0,0 

POINT 186 O.IOOOE Ol 
0,0 
0,0 

POINT 185 O.IOOOE Ot 
0,0 
0,0 

POINT 178 O.IOOOE 01 
0.0 
0,0 

POINT 170 O.IOOOE 01 
0,0 
0,0 

POINT 163 O.IOOOE 01 
0,0 
0,0 

POINT 171 O.IOOOE 01 
0,0 
0,0 

POINT 166 O.IOOOE 01 
0,0 
0,0 

POINT 187 O.IOOOE 01 
0,0 
0,0 

POINT 199 O.IOOOE 01 
0,0 
0,0 

POINT 168 0,1 OOOE 01 
0,0 
0,0 

POINT 163 O.IOOOE 01 


0, 

0 


0.0 


0* 

1 OOOE 

01 

0,0 


0, 

0 


O.IOOOE 

01 

0. 

0 


0.0 


0* 

1 OOOE 

01 

0.0 


0* 

0 


O.IOOOE 

01 

0. 

0 


0,0 


0, 

1 OOOE 

01 

0,0 


0. 

0 


O.IOOOE 

01 

0 . 

0 


0,0 


0. 

1 OOOE 

01 

0.0 


o. 

0 


O.IOOOE 

01 

0. 

0 


0,0 


0. 

1 OOOE 

01 

0.0 


0, 

0 


O.IOOOE 

01 

0. 

0 


0,0 


0. 

1 OOOE 

01 

0.0 


0, 

0 


O.iOOOE 

01 

0, 

0 


0,0 


0* 

1 OOOE 

01 

0,0 


0* 

0 


O.IOOOE 

01 

0. 

0 


0.0 


0* 

1 OOOE 

01 

0.0 


0* 

0 


O.IOOOE 

01 

0, 

0 


0.0 


0* 

1000E 

01 

0.0 


o. 

0 


O.IOOOE 

01 

o. 

0 


0.0 


0* 

1 OOOE 

01 

0.0 


0* 

0 


O.IOOOE 

01 

0. 

0 


0.0 


0» 

1 OOOE 

01 

0.0 


0* 

0 


O.IOOOE 

01 

0, 

0 


0.0 


o. 

1 OOOE 

01 

0.0 


0. 

0 


O.IOOOE 

01 

0. 

0 


0.0 


0. 

1 OOOE 

01 

0.0 


0. 

0 


O.IOOOE 

01 

0, 

0 


0,0 


0. 

1 OOOE 

01 

0,0 


0. 

0 


O.IOOOE 

01 

0. 

0 


0.0 


0, 

1 OOOE 

01 

0,0 


0, 

0 


O.iOOOE 

01 

0* 

0 


0,0 
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POINT 217 


0,0 

0,0 

O.lOOOE 01 
0,0 
0,0 


O.ICOOE 01 

0.0 

0.0 

o.i root os 
o.o 


o.o 

0.1000E 01 

0,0 

0,0 

O.lOOOE 01 
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I I F PAT ION NUMBER 


1 



* SIGHA ZERO* 0.1186E-01 * 



transformation pApA|METtR$ 


SCALE FACTOR 
translation IN 
TRANSLATION IN 
TRANSLATION IN 
ANGLE OMEGA 
ANGLE PHI 
ANGLF CApA 


OLD VALUE 

0.CP667 
X 10088.594 

Y S9«?.703 

2 434.501 

0 0 46,25 

0 5 41 ,50 

0 55 29,80 


CORRECTION 

-0, 1243e-05 
-0,133*E 04 
0,«ftri5F 03 
-0.2230E 02 
0 “2*23,25 
0 “3-27,10 
0 1 6,19 


NEW VALLE 

0,00667 
8754,570 
6463.172 
412,202 
0 -1-37,00 
0 2 14,40 

0 56 35.99 


ROTATIONS IN RAnlANS 


ANGLE omega 
angle PHI 
ANGLE CApA 


0.00022424 
0.0Q lf 5564 
0,01614332 


-0.6945p-03 
-0*1 OC 4^-02 
0.3209F-03 


-0*00047026 

0.00065161 

0.01646423 


DEVIATION 

0,2846e-05 
U.7299E 01 
0.AB96E 01 
0.1213E 02 
0 3 52.87 
P 8 46,98 
P 1 54,37 


P. 1 129 e -o2 
p,2555f-o? 
p ,5545 f-03 



the solution converges after iteration number < 



SIGMA ZEROS 0 . 1187 E -01 * 



TRANSFORMATION PARAMETER*; 


OLD VALUE CORRECTION FINAL VALUE 


deviation 


SCALE FACTOR 
TRANSLATION IN x 
TRANSLATION IN Y 
TRANSl ATION IN 7 . 
ANGLE OMEGA 
ANGLE PHI 
ANGLE CAPA 


0*006*7 
8754, 57Q 
6463.172 
4 1 ? o ?02 
0 -1-37.00 
0 214.40 

0 56 35*99 


0 , 10 * 7 F-U 7 
- 0 . 407 3 f 00 
- 0 . 3552 F 00 
• 0 . 1665 F 01 
0 0 0.66 
0 0 1.55 

0 0 -O.O 5 


0.00667 
8754,160 
6462,816 
410,537 
0 - 1 - 3 6 ,12 
0 2 15.95 

0 56 35.94 


O . 28 /) 7 E "0 5 
0 ♦ 7 3 3 ? E 01 
n. 4957 f oi 

P, 1?59f 0? 

0 3 53 , 1 a 

0 8 47,4? 

O 1 54 . 4 1 


THE VApI ANCE-CDVA R lANCE MApTlX OF THE TRANSFORMATION 


0 , 81 o 5 p. 

e. i 6 & 3 f c 

0 « 4 4 ?9 p" 
0 * 8866 *. 
0 * 6864 *- 

0 * ? 1 06p* 
0 * 1 069 p- 


11 0 , 1683 *. 

05 0.5376* 

06 0 . 1 055 p 
07 - 0 , 1 37 2 E 
10 0 . 2339 F" 
1 0 * 0 . 1 1 47 p* 
09 - 0 . 169 ]*- 


05 - 0 , 4428* -06 

02 0 * 1 055 p on* 
00 0*2457 p 0 ? 
00 0 * 8379 f *oi 
04 0» 8 1 6 1 f* 05* 

03 O* 6068 f*O 4 
03 * 0 * 4052 F- 03 ' 


0 , 88 7 op« 
0 . 1 37 2 (r 
0.838 If* 
Ot 1584 c 
0 . ?888 *■ 
0 * 84 55 p' 
0.6806F* 


07 0 . 6864 E -10 
00 0 . 2 339 f- o 4 * 
0 ] 0 » 8 1 6 1 £- o 5 
03 - 0 . 2887 E -03 

03 0 . 1278 e- 05 ' 
02 - 0 . 6755^-06 

04 0.1128e-o7« 


0 , 2105 *- 
0 * 1 1 47 l* 
0* 6086*' 
0 . 6453 L- 
0 * 67 55 E* 
0 * 6538 ** 
0 . 5468 ** 


10-0 , IO 69 F-O 9 

03- 0* 169 l f" 03 

04 - 0 . 4052 E ”03 
' 02 - 0 * 68 o 4 f -04 

06 0 * 1 1 ?6 F“ 07 
05"0*5468f"07 

07 o» 3077**06 
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LIST OF ADJUSTED GROUND COORDINATES 


NUMBER 

X 

Y 

Z 

SIG(X) 

SIGCY) 


SI6C Z) 


27 

1«0?.769 

15105, 105 

29 | , 357 

0, 9756 E 01 

0, 76«6E 

01 

0.2450F 

0? 

?9 

10746,65? 

21621, 32c 

392,987 

0,11 35f 02 

0 . 8 3 1 4 r 

01 

0 , 2 1 82F 

02 

IS 

14106, 6 29 

-9589,346 

8 8 4 , 3 1 8 

0«ll57r 02 

0.6437E 

01 

0.2166E 

02 

1 4 

6000,916 

•6094, 60S 

384,560 

0* 99 73 E 01 

0.7317E 

01 

0, 1992E 

02 

17 

8703,076 

5197,234 

305,567 

0.7485F 01 

0*51 75f 

ol 

0. 1281F 

02 

30 

1722?. 941 

9636.289 

349,100 

0* 6229 p 01 

0. 66«4p 

01 

0 , 23q4 F 

02 



LIST 

Of COMPUTED 

GHCl'NO COORDINATES df 

UNKNOWN POINTS 




NUMBER 

y 

Y 

7 

SIGMA X 

SIGmA y 

SIGMA 2 

30 

1 7223 ,25g 

9636,402 

350.913 

p.466 

7.122 

25,347 

7 

16100. 777 

12002.691 

928,696 

0.731 

7.025 

23.511 

? 8 

2 7 3l .TO 7 

2 2347.523 

345.682 

1 1 * 9 ?2 

9.163 

26.839 

ITS 

7848.547 

-11888.434 

462.655 

12.623 

9,404 

24.422 

1 8 ? 

ll9l6.59o 

'11058, H3 

882.592 

12.438 


24.868 

190 

17537,324 

-10721 .3*5 

1048.351 

12.903 

10.241 

32,288 

l»t 

11899.895 

-6171.297 

709.674 

10.403 

7,605 

20.721 

174 

7330.199 

-6822.988 

366.749 

10.500 

7.72} 

19.988 

166 

2593. 102 

'6628,355 

388,037 

10.676 

0.351 

25.058 

172 

7165. 727 

3173. 4f 6 

327.369 

7,729 

5.439 

13.820 

1T3 

69o7. o«7 

-2153. 035 

362.430 

8.896 

6. 4 34 

16.669 

180 

11896.601 

-1012.883 

422 . 066 

8,687 

6.320 

17,224 

1 88 

17463.773 

- 1767,74? 

737,714 

9.633 

7,914 

27,262 

169 

T128.383 

17310.672 

362. 67t 

9,644 

7,005 

18.120 

ITT 

1 1690.988 

17778.391 

385.775 

9,811 

7.267 

19,511 

1 8 4 

16753.371 

10040 . 121 

868,005 

10,411 

8.437 

27.392 

185 

17417.273 

13783,797 

840. 036 

9.279 

7,740 

26.830 

ITS 

12537.379 

14335.777 

410.233 

8.778 

6.521 

18,260 

1 T 0 

7025.543 

14157.277 

374.324 

8 , 66 i 

6.207 

15.996 

163 

2497.0*6 

7 115,887 

314,164 

7.964 

6 . 242 

20.426 

l 7 1 

79 ? 2. l 7p 

9?47.543 

392.542 

7,653 

5.3 34 

13.233 

186 

]69ftp.68o 

8644,910 

372,404 

0.328 

6 • 95 1 

24.599 

1 «T 

16739.922 

4478. 9i4 

483,951 

0.319 

6.860 

24.138 

199 

13157.523 

2976.610 

233. 7 fl ? 

7.988 

5.920 

17.406 

168 

68o4. 008 

22749, 4*7 

399,666 

1 1 ,784 

8.736 

22.903 

183 

17216. 0®6 

23514.016 

786,261 

12.508 

10.100 

31.629 

21 7 

12064.441 

23087.430 

453.056 

11 • 9 1 0 

8.975 

24.197 
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6,2 Tn STUDY accuracy of absulute orientation by SIMULATION 


6,2,1 INPUT OATA 


1 

,01 

25 

25 e 

• OOOF +06 

1.0 

1 , 000 F + 06 1 . 000 F + U6 1 . OOOF + 06 
1 

- 1000,0 

1000.0 

745,0 

- 500, 0 

1000,0 

730,0 

.0 

1000,0 

747.0 

500,0 

1000.0 

735,0 

1000,0 

1000,0 

780,0 

- 1000,0 

500,0 

790,0 

- 500.0 

500,0 

800,0 

.0 

500,0 

720,0 

500.0 

500.0 

735,0 

1000.0 

500.0 

749,0 

■=■ 1000,0 

.0 

713,0 

- 500,0 

.0 

784,0 

.0 

♦ 0 

795,0 

500,0 

.0 

810,0 

1000,0 

.0 

725.0 

- 1000.0 

- 500.0 

730,0 

- 500.0 

- 500.0 

747,0 

.0 

- 500.0 

758.0 

500.0 

- 500.0 

739,0 

1000,0 

- 500,0 

729,0 

- 1000.0 

- 1000,0 

786,0 

- 500.0 

- 1000.0 

793,0 

.0 

- 1000.0 

749,0 

500.0 

- 1000,0 

763,0 

1000,0 

- 1000,0 

777,0 

1 » 44 E +02 

1 , 44 E +02 

25 E + 0 ? 

1 . OOE+OO 

l . OOE+OO 5 . 00 E +00 
200.0 - 200.0 25.0 

43211 

200.0 

1431 15 

1?3455 

479213 

69813 1111 

.too 

? 0*0 -15 

.0 5.0 


1 » OOOF + 06 1 * OOOE+06 1, 000F+06 1. OOOF+06 


•33142 


-.33142 ,67189 


12356791 213456789 

.01 -0.0? -0,03 
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6.2.2 OUTPUT 

GENERAL DATA 

NUMBER OF POINTS USEE IN THE SOLUTION 25 

MAXIMUM NUMBER OF ITERATIONS All OWED 8 


THF VARI ANCE-CnVARlANCF MATRIX FOR THE TRANSFORMATION PARAMETERS 


(ORDER OF VARIABLES * SC ALF* DX* PY, OZ > OME> PHI » C APA ) 


0 , 1 OOF 

07 0,0 

0,0 

0,0 

0,0 

0,0 

0.0 

0,0 

0,1 OOF 

07 0,0 

0,0 

0,0 

0,0 

0,0 

0.0 

0,0 

O.lOOE 07 

0,0 

0,0 

0.0 

0,0 

0,0 

0,0 

0.0 

0,1001 

07 0,0 

0,0 

0*0 

0,0 

0.0 

0.0 

0,0 

0 . 100 E 

07 0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0 . 100 E 

07 0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0 . 100 F 07 



* SIGMA ZEROa O.iOOOO * 



83 

*********************** PRINT out or thf 


Ll$T or TRU * coordinates 

POINT 

NUMBER 

Ur POINTS USED 
X 

in the solution 
coordinates 

Y 


7 

t 

•■ 1000,000 

1000,000 

745 

*000 

2 

- 500,000 

1000,000 

730 

.000 

3 

0,0 

1000,000 

747 

,000 

4 

500,000 

1000,000 

735 

,000 

5 

1000,000 

1000,000 

780 

,000 

4 

- 1000,000 

500,000 

790 

.000 

7 

- 500,000 

500,000 

800 

,000 

a 

0,0 

500,000 

7?0 

,000 

9 

500,000 

500,000 

735 

.000 

10 

1000,000 

500,000 

749 

,000 

11 

- 1000,000 

0,0 

713 

,000 

12 

- 500,000 

0,0 

784 

,000 

13 

0,0 

0,0 

795 

.000 

14 

500,000 

0.0 

810 

*000 

15 

1000,000 

0,0 

725 

.000 

16 

' 1000,000 

- 500,000 

730 

*000 

17 

- 500,000 

- 500.000 

747 

.000 

18 

0,0 

- 500.000 

758 

,000 

19 

500*000 

• 500,000 

739 

.000 

20 

1000,000 

- 500,000 

7?9 

.000 

21 

- 1000,000 

- 1000,000 

786 

.000 

22 

- 500,000 

- 1000,000 

793 

.000 

23 

0.0 

- 1000,000 

749 

,000 

24 

500.000 

- 1000,000 

763 

.000 

25 

1000,000 

• 1000,000 

777 

,000 
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THF VARlANCfc-COVARl 


VARIANCF-CnvARtANCE MATRICES ******************** 
ANCE MATRICES FOR MODEL COORDINATES 


0,1«AE 03 

0.0 

0.0 


0.0 

O.lAAt 03 

0.0 


0.0 

0.0 

0,2?0E 02 


THE VARlANcf-COVARIANCE mATPIcES FOR GPOUNl) COORDINATES 


O.lOOf 

0.0 

0.0 

ThE ratio sIgm a 
t h e ratio SIgm a 
t h e ratio sIgm a 


1 0.0 

0.100E 01 

0.0 

GROUND TO Sltih*A 
GROUND TO sIgm a 

ground Tn sIg^a 


0.0 

0.0 

0.500F 01 
MODEL, e OP TrE X 
mOqEl, top TmE y 

model, f or t h e z 


COORDINATES, at 
COORDINATES, at 
COORDINATES, Af 


the ground sc a l* 
the ground sc A Lf 

ThE GROUND Sc a L* 
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OUTPUT or SUBROUTINE INPUT1 


transformation parameters 


PARAMETER QENf rated corrections approximations 

SCALE 1,00000 0,10000 1,10000 
TRANSLATION IN X *00.00000 20.0000 220.0000 
TRANSLATION IN Y "200,00000 "15.0000 -215.0000 
TRANSLATION IN Z 25,00000 5.0000 30.0000 
ANGLE OME 0,331*2 0,01000 0.341*2 
ANGIE PHI -0,33142 -0,02000 -0,3514? 
ANGLE CApA 0.67189 -0,03000 0,6*189 

ANGLES IN DEGREES 

ANGLF ONEGA 18 59 20.23 0 34 22,65 19 33 42.91 
ANGLE PHI -18-59-20.23 -1 -8-45.29 -20 -8 -5.54 
ANGLE CAPA 38 29 47,20 -1-43 -7,94 36 46 39' P6 


MODEL COORDINATES 


POINT 



GENERATED 




pERl URpEn 



number 

X 


Y 

7 


X 

Y 


Z 


1 

22647. 

64 

10898,97 

55011 , 

57 

22660.36 

18926, 

43 

55018 

34 

2 

27240. 

86 

13060.09 

40730, 

49 

27228,74 

1306'» 

63 

48726 

67 

3 

32529, 

18 

1 lV?4,55 

51076. 

73 

32520.73 

11929. 

30 

51064 

63 

4 

33359. 

00 

8150.96 

4765°. 

56 

33369,06 

816 ?. 

49 

47658 

86 

5 

39J74. 

9 1 

5230. 15 

50240. 

33 

39279,86 

5238. 

57 

50237 

13 

6 

19429, 

45 

11586,27 

5P 169. 

51 

19417.36 

11573, 

68 

58 1 6 8 

65 

7 

30584. 

44 

1 1 9 1 5 » n9 

56303. 

95 

30586.27 

11915. 

05 

5638 9 

55 

8 

26841* 

99 

6459, 16 

5 1 7 4? , 

48 

26862.22 

644?. 

11 

51747 

3? 

9 

32007. 

84 

2993.38 

46519, 

10 

32000.60 

3016* 

47 

48513 

02 

10 

36673. 

04 

96P. ’1 

4ftl49, 

16 

36662.23 

9rf. 

30 

48146 

82 

11 

17534, 

02 

1 19?9,2o 

52396, 

84 

17555,44 

11931. 

34 

52392 

10 

12 

240* 9 * 

79 

8177 . 08 

57556. 

79 

2408 1. 67 

8i9r. 

9 1 

57562 

21 

13 

29 184 . 

73 

5212.50 

5657 r . 

44 

?9 164,67 

5203. 

15 

56576 

87 

14 

32848. 

O 6 

2 9 3 9 « 8 l 

54576, 

1 0 

32839,83 

2947, 

67 

54575 

05 

15 

30021* 

73 

-2274 , 08 

50215. 

87 

30023, 14 

•2?7 * • 

26 

5 0207 

72 

1« 

15576. 

11 

4449,33 

5»39o. 

55 

15562.32 

4461 , 

25 

58401 

98 

17 

20033. 

22 

5841.73 

55019. 

64 

20030.57 

5850* 

76 

55022 

61 

1» 

24567. 

80 

6056 .96 

54833. 

3? 

24553.76 

6026 . 

88 

54035 

97 

1 9 

29265. 

7ft 

-2974,55 

50120. 

10 

29244.99 

*2966, 

24 

501 12 

75 

20 

33212, 

50 

-522P .58 

49567. 

02 

33219,15 

-5237, 

89 

49565 

62 

21 

19868. 

9 1 

1 3?? , 58 

6 0887 , 

37 

19873.79 

1315. 

47 

60884 

85 

22 

19702. 

75 

1277,97 

60456. 

94 

19691, o7 

1275, 

43 

6045 9 

58 

?3 

23498, 

14 

-4775 , 02 

54550. 

66 

23400* 54 

- 47p f . 

35 

54561 

93 

24 

25548. 

37 

-2776 » 49 

58456. 

63 

25543.25 

-2766. 

10 

58454 

64 

25 

31565, 

46 

-12087, fll 

54010. 

04 

31542.93 

-1210*. 

13 

5400 9 

25 
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1 O 
I X 

a 

UNO CO 

ORPIN 

A T E S 



POINT 


genfrated 



RFRU'RBFp 


NUMflFR 

X 

Y 

7 

X 

f 

Z 

1 

•12063.27 

9559,49 

60446,23 

-12064,34 

9559.05 

60448.75 

7 

•3183.52 

9io9, 07 

56497.87 

•3183, o7 

9110.84 

56499 , 38 

3 

635.Q6 

10146.43 

60868,55 

635,86 

10148.25 

6o«67.59 

a 

#562. 18 

8576.1 3 

57942.81 

4583,43 

6574,65 

57945.0? 

5 

9839.71 

8419,86 

62688,69 

9839,74 

8419,20 

6?6P9. 0 ? 

6 

"11168.25 

1066,59 

61381,31 

-11166.48 

1065,62 

61380.26 

7 

-?5?5.48 

7522. 38 

64751.07 

-2525.70 

7521.35 

64747*08 

8 

-573.51 

2660.53 

58597 . 32 

•573,62 

2661.73 

58597,84 

9 

6338. ?9 

3#71.53 

57786,64 

6339. 35 

3470.84 

57784 , 39 

m 

11 102.75 

4313.04 

59396, 3 ^ 

11101.42 

4312.84 

§9398.19 

it 

-10695.13 

2159,91 

55402.93 

•10894,72 

2159,19 

55401 *74 

12 

-5528.15 

854 • 1 8 

62675,47 

-5527 » 18 

85«.26 

6?674 » 9 l 

13 

323. oT 

1355*10 

63874.04 

323.47 

1356.92 

63870. 1 1 

1* 

5 020* 1 1 

j989 , 06 

63573,26 

5021. 16 

1989,59 

63573.28 

15 

7416,74 

-2300*33 

58o9fi,63 

T4 18 » 1 1 

•2301*87 

5«o91.58 

16 

-989 i , 12 

-67o?.27 

59419,94 

-9890,89 

-6701*67 

59422.27 

IT 

-6315.31 

"2288.80 

58469,88 

-6315,28 

*2290*86 

58471. 04 

1 ® 

-3o?5,53 

235,26 

6o32T.?7 

-3025,05 

2*5.55 

6o3?4 , 16 

l 9 

7300* 75 

- 32 1 7 , 7 2 

57625,63 

7301.79 

-3214.8? 

576?3.96 

?ft 

11728.21 

-26t8.?4 

58738,54 

1 1728.29 

"2866.50 

58737,78 

21 

-5686.25 

-7Bj9,5 3 

63358 . 16 

-5684.55 

• 78 1 9 , 1 1 

63361.55 

?? 

"5643s 56 

- 7 8 q7 « 7 q 

62698.6, 

-5643. 34 

* 78 06 , 8 0 

6?899. 0 9 

23 

?647,53 

*8951 , 65 

589j4 • 4 1 

2649,oi 

-8950*93 

58927, 0 4 

?# 

1 ^ 20 » 48 

-75 r 3.o7 

63443.93 

1 7 1 9 . 1 1 

-7500*95 

63444.42 

25 

13100» 9 0 

- l q596 , $9 

61547 , 85 

13101*68 

-J0595.79 

6i547.82 
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WEIGHTS USED IN SOLUTION 


ft********************* 


TmF w EIqhT matrix for The TPANSfOPmATIPN PARAMETERS 


0 , 1 O 0 F -07 

0.0 

0,0 

0,0 

o.n 

0,0 

0.0 

0.0 

O . I 00 E - O 7 

0,0 

0,0 

0,0 

0.0 

0,0 

0.0 

0,0 

o . iooe’ot 

0,0 

0,0 

0,0 

0,0 

0.0 

0.0 

0,0 

0 . 100 E -07 

0,0 

0,0 

0,0 

0.0 

0,0 

0,0 

0,0 

0 . 100 E -07 

0.0 

0.0 

0.0 

0,0 

0,0 

0,0 

0.0 

o . ioor -07 

0,0 

0.0 

0,0 

0,0 

0,0 

0,0 

0.0 

0 . I 


thf ground coordinate weight matrix 


O . lOOOf *01 
0,0 
0,0 


0,0 

0 , lOOOF-OI 

0,0 


0,0 

0.0 

0 , 2000 E *02 


THE MODEL COORDINATE WEIGHT MATRIX 


0.69MF-04 

0,0 

0,0 


0,0 

0, 69MF*04 

0.0 


0.0 

0,0 

0 . 4000 E -03 
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JURATION NUMBER 


1 ft*********************** 


+++++++********+*****+«++******* 
* * 


♦ SIGMA ZERO* O f i0^3E 00 * 



TRANSFORMATION PARAMETERS 


OLD VALUE CORRECTION 


NFN VALUE DEVIATION 


SCALE F A C 1 OR 




1,10000 

- 0 , 1000 E 00 

0.9999b 

translation 

IN 

X 


220.000 

0.3248E 02 

252,477 

translation 

IN 

Y 


-215,000 

-0.2085F 02 

-235.852 

translation 

IN 

7 


30,000 

0 , 98 A 3f 01 

39,883 

ANGLE ONEGA 



19 

33 42.91 

0-31-20.22 

19 2 2?,66 

ANGLE PHI 



-20 

-8 -5,54 

1 1 47,42 

-19 -6-18,09 

angle CApA 



36 

46 39*26 

1 35 28,29 

36 22 7,53 


C.2122E-03 
0.1065E 0? 
0 • 1 2 2 3 E 0? 
0*ll75E 02 
0 0 39.80 
0 0 31.35 
0 0 52 , ?4 


POTATIONS IN RAn I AN5 


angle OmEGA 

angle phi 
angle CApA 


0,34141994 

-0.35U1993 

0,64188993 


- 0 . 9116^-02 

0»1797e-01 

0.2777e*01 


0.33230430 

-0.33344579 

0.66966146 


0 . 1934^-03 

0.1520E-03 

0.2533E-03 
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I I f FAT I DN NUMBER 2 


T P A N 


SCALE FACTOR 
translation 
translation 
translation 
angle omega 

ANGLE PM I 
ANGLE CAPA 


POTATIONS in 


angle OMEGA 

ANGLE PHI 
ANGLE CAPA 


******************************* 

* 

* SIGMA ZERO* 0. 1031E 00 

* 


SfORMATTON PARAMETERS 


OLO VALUE CORRECTION NEH VALLE 


DFvIATIPN 


0.99998 
IN X 252.47 7 

IN Y -2^5.852 

IN Z 39.883 

19 2 22,66 

-19 -6“1 P • 09 
38 22 7,53 


0,55f 1F-03 
•0.33fl9f 02 
0.2131E 02 
0.1961F 02 
0 -3-26.fi 
0 6 34.05 

0 9 11,26 


1 ,00053 
218.986 
-214.538 
59,494 
10 58 55.95 
-18-59-44.07 
38 31 18,83 


P.2098E-03 
P.1134E 0? 

O. 1300E 02 

P. 1278E 02 

P 0 42.41 
0 0 33.41 

f 0 56,42 


RADIANS 


0. 33230430 
“0,3334*579 
0,66966146 


-0. 1 002f"02 

0. 191 OF-02 
0. 2673p-02 


0.33130211 

-0.33153540 

0.67233402 


0.?056e-03 

O. 1620E-03 

P. 2735F-03 
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FINAL OUTPt'T 


ThF SOLUTION CONVERGES AFTER ITERATION NUMBER 3 



* SIGMA ZERO* 0 * 1 03 1 E 00 * 


* * 


TRANSFORMATION PARAMETER*; 


OLD VALUE CORRECTION FINAL VALUE 


deviation 


SCALE FACTOR 
TRANSLATION IN 
TRANSLATION IN 
TRANSLATION in 

angle omega 

ANGIE PHI 
ANGLF CAPA 


1 ,00053 
X 218.986 

y -214.518 

l 59.494 

18 58 5?. 95 
-16-59-AA.07 
38 31 1«.83 


0 ,54<S4e»05 
-0.3506F 00 
0.1443F 00 
0. 1 574 £ 00 
0 0 0»02 
0 0 0.01 
0 0 - 0.46 


1 ,00054 
218.635 
-214.393 
59,651 
18 58 55.95 
-18-59-44,01 
38 31 18.33 


0,20 9 9E-03 
0.1I32F 02 
0.1297f 02 
0,J?77j; 02 
C 0 42,34 
0 0 33,37 
0 0 56,36 


the VArIANcE-cOvArIAMcE mArTIx Of the TRANSFORMATION 


Oo 44 0 4 f -o7-Os 6 5!2£-o 3-Q.6 344 f - 0 3 q, ?66 j f * 02* Q» 1?10E"Q 7 0» 1 q33e - 0 7 *0» 335 1 j- 08 
-p.65i 2f-o3 0.128 of 03“O»3355f 00*0. 3776^ 02 0. 5656£- o3*0* 1 58 0£* 02 0»17 o7f- 0? 
-O.8344 f -O3*0. 3352f 00 0* 1*6 3p 03-0.38f5f 02 0. 2380E* 02" 0. 339 ?£- 03- 0. 9976 z ” 03 
0» 266 1 f* 02- 0. 3776 f o 2 ” 0 * 3865 f 0 ? 0 • 1 6 3 1 E 03- 0. 7557 E - o3 0. 62*U W 03"0* l9o6f-03 
-0»1210F”07 0.5656^-03 0» 238 0 f*02- 0* 7557 £-03 0* 4214e-07-O. 7335e"08 0»T780£-08 
0. 1033f-07-0. 158 of* 02-0. 3391 F -03 0*6 26 2e-03-0. 7 334e*o8 0. ?6 17e*’07-o.4897f-08 
- 0 » 335 1 F" 08 0, 17o7f*02"0.9976f-03-0»1906f-03 0. 778 OE" 08- 0. 4697 E* 08 o»7465e"07 



LIST OF ADJUSTED GROUND COORDINATES 


NUMBER 

X 

Y 

Z 


5 1 G ( 

K ) 

STG(Y) 


S1G(Z) 


1 

•12064,352 

9559,166 

60449,133 

0, 

1027f 

01 

0 • 1 027 E 

01 

0 » 2 1 7 3 E 

01 

2 

•3163, 1©7 

9110.613 

56498,746 

0, 

!026f 

01 

0* 10?6E 

01 

0.2158F 

01 

3 

635, 91B 

10145.344 

60865, 816 

0. 

1026e 

01 

0. 1 0?6E 

ot 

0.2159E 

01 

4 

4583, 4 0 6 

8574.953 

57945.219 

0* 

1026E 

01 

0* 10?6f 

01 

0 , 2 1 56 E 

Ol 

5 

9839,723 

8419,285 

62688.859 

0, 

1026E 

01 

0.1C27E 

01 

0 » 2 1 7 1 E 

01 

6 

•11166.414 

1065.517 

61379,6*6 

0* 

1026E 

01 

0*1 026f 

01 

0.2157E 

01 

7 

-2525,75? 

7521,273 

64748,109 

0, 

1026E 

01 

0. 1026F 

01 

0.2159E 

01 

6 

-573.482 

266 1,634 

58598.977 

0* 

1026E 

01 

t* 1 U?6 E 

01 

0» 2 1 43 E 

Ol 

9 

6339.230 

3470,96? 

57784.301 

0* 

1026E 

01 

0. 1 026E 

01 

0.2149E 

01 

10 

11101,340 

4312, 8*0 

59397,992 

0, 

1 02 6E 

01 

0.102 6 E 

01 

0 , 21 5 9E 

01 

11 

-1 P 894,473 

2159,314 

55401 ,160 

0. 

1026E 

01 

t * 1 0?7 E 

01 

0 * 2 1 67 E 

01 

1? 

-5527, 2o7 

854.36® 

6j>675»7?7 

O' 

1 026f 

01 

0»10?6e 

01 

0.2147E 

01 

13 

323,315 

1356,696 

63871,277 

0, 

1026E 

01 

C»l0?6r 

01 

0.2146E 

01 

14 

5021.102 

1989,643 

63572.836 

0* 

1026e 

01 

0»10?6e 

ol 

0*2l5oE 

01 

15 

7418,184 

•2301 , 76o 

56o90,805 

0, 

1026E 

01 

0» 1 0?6 e 

01 

0.2150E 

01 

16 

-9891,016 

-67oi*73o 

59423.3M 

0» 

1026E 

01 

0» 1 0?7 e 

01 

p. 21 66E 

01 

17 

-6315,261 

-229o,8 16 

58471,195 

0' 

1 026e 

01 

0*1 026f 

01 

0.2151E 

Ol 

1® 

-3025.025 

235.296 

6o324,7*0 

0. 

1026e 

01 

O' 1 0?6f 

01 

0.2143E 

01 

19 

73ol .754 

-3214.919 

57623,543 

0* 

1 026p 

01 

0* 10?6e 

01 

0»2152E 

01 

20 

11728.320 

-2866.526 

58738.539 

0. 

1026E 

01 

fc» 1 0?6f 

01 

0.2160E 

01 

21 

-5684,367 

-7819,023 

63360.463 

0* 

1 026 E 

01 

0. 10?6E 

01 

0.2159E 

01 

22 

-5643,301 

•7806.816 

62898,664 

O' 

1026E 

01 

0 • 1 0?6 E 

01 

0 • 2 1 56 E 

01 

23 

2648,951 

•8951,008 

58927,180 

0' 

1026E 

01 

0* 1 026f 

01 

0.2159E 

01 

24 

1719,12? 

-7500.844 

63443,965 

0, 

1026f 

01 

0* 10?6e 

01 

0 • 2 1 54 E 

Ol 

25 

13101.598, 

“10595.902 

6i548»op8 

0' 

1027E 

01 

0*1 027f 

01 

0* 218 1 E 

01 
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6,3 TO DETERMINE THE UNCERTAINTY IN THE ORIENTATION OF THE SURFACE DEFINED 


PY A SET DF TRIANGULATED PASS POINTS 


6.3,1 INPUT DATA 


3 

1.0 

9 

9 8 

OOOE+06 

1, 000E+06 1 

♦000F+06 1. OOOE+06 

12345 

1234567 

123 54321 

7 

•205573.4 

•’84998,6 

8 

-183918.4 

•8944} ,6 

9 

-161431*8 

-94340.5 

49 

-227384.6 

•27526.5 

50 

-204120. 1 

-32641.7 

51 

-179883.2 

-38225.9 

91 

-250402*5 

35956.8 

92 

-225748.3 

30354,6 

93 

-199706.4 

24166.2 

7 

6.4 

11.3 

8 

7,1 

19,7 

9 

6.3 

10,3 

49 

6,4 

5,3 

50 

8,0 

6,6 

51 

6.3 

5,4 

91 

6,6 

10,4 

92 

11*4 

15,6 

93 

7,6 

9,3 

0.5 

0*1 


50.0 

10,0 


0,5 

0,2 


1.0 

0.2 


OO0EM 0 

1 ,000E 

-10 1.000E-10 


1 1 

.OOOE+06 1 , 000E+06 1. OOOE+06 

84675.1 

86030.0 

89462.9 
84682*8 
86455.7 
9 1 48 7t 4 

80707.9 
82398.3 
86004*7 

16.9 

31.2 

14.6 

16.4 

29.4 

13.5 

20.3 

31.7 

17.7 9S 
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6.3.2 OUTPUT 

GFMtRAL OATA 


9 

8 


NUMBER OF POINTS USfP IN THE SOI UT fON 
MAXIMUM NUMBFR OF ITERATIONS All OWFU 


THF vARlANCf *C0VARIANCF M A TPly FOR THE! TRANSFORMATION PA&AmUFRS 


(•ORDER OF VARIABLES = $ C AL F * OX# P Y» f)Z, QMF » PHI t C AP A ) 


0.100F 

07 0,0 

0,0 

0.0 

0,0 

0,1 00e 

07 o.o 

0,0 

0.0 

0.0 

0 . 100 E 

o 

• 

o 

© 

0.0 

0,0 

0.0 

0 , 1 OOF 

0,0 

r.o 

0.0 

0,0 

0,0 

0,0 

0,0 

0.0 

0.0 

0.0 

0.0 

0,0 


0,0 0,0 0.0 

0,0 0.0 0.0 

0,0 0,0 0.0 

07 0,0 0,0 0.0 

0 , 1 OOE 07 0.0 0.0 

0,0 0 , 1 OOF 07 0,0 

0,0 0,0 0,1 OOf 07 


* 

* 


* 

SIiiM ZE . RO * 1,00000 * 
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THF GTVFn MODEL COURDInMES OF THE CONTROL POINTS 


DUMBER 

1 7 

? ft 

3 9 

4 49 

5 50 

6 51 

7 91 

e 9 ? 

9 93 


X 

" 205573.39 
-183910.39 
-161431.75 
- 227304.56 
-2041^0.06 
-179083.19 
-250402,50 
-225746. 25 
-199706.39 


Y 

-0/i998 » 56 
-69441 .75 
•94340.50 
-27526*50 
-32641.70 
-30220*90 
3*950.80 
30354.60 
24100*20 


7 

64675*06 

66030.00 

69462,80 

64662.75 

66455.69 
91467,38 
b07L7,88 
1 23^9 . 25 

66000.69 



O ® N O U1 frw f\) 
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THF DEVIATIONS 

OF THE MOpEt 

COORDINATES UE 

THE CONTROL points 

NUMBER 

SIGMA X 

SIGMA Y 

SIGMA Z 

7 

6,3999996 

1 1,2909992 

16.8999939 

s 

7,0999994 

19,6909969 

31.1999969 

9 

6 , 2999992 

10.2909992 

14,5999994 

A 9 

6,3999996 

5,2999992 

16,3999939 

50 

8,0000000 

6.5909994 

29.3999939 

51 

6,2999992 

5.3909906 

13,5000000 

91 

8,5999994 

10,3909996 

20,299987b 

9? 

11 , 3999996 

15.5909904 

31,6999969 

93 

7,5999994 

9.2909902 

17.6999969 
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*********************** PRINT OUT UF !NPI'T3 


THE ARBITRARY ROTATIONS OF THE MODEL 


IN R A U IN DEGREES 


DmEGA 

o ,0067650 

0 

23 

16,171 

PHI 

0,0116747 

0 

AO 

A9.331 

CAPA 

0,0116670 

0 

AO 

A 7 ,901 


APPROXIMATIONS FOR T HE 1 RANSfORMAT If|N pARAMFTtRS 


PARAMETER 


APPROXIMATION 

in pad in degrees 


SCALE 

translation IN X 
translation IN Y 
translation in z 

ROTATION OMEGA 

rotation phi 
POTATION CApA 


1 .04620 
39, 48483 
47, 34769 
27,77567 

-0,00326 0 -11 -12,320 
*0,01163 0 *39 -59,?50 
-0,01065 0 -36 -36,780 
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ground coordinates * (gyrated by rotation of the given mooel coordinates) 



NUMBER 

X 

Y 

Z 

1 

7 

-207556.49 

-81965.94 

82801.31 

2 

e 

*185975,50 

-86656,38 

84443. 38 

3 

9 

-163591.38 

-91797.56 

68196.06 

A 

49 

•228680.00 

-24240.52 

62161 .06 

5 

So 

-205500.81 

•20619.05 

64240.69 

6 

51 

-181393,38 

-35455.57 

69601 .50 

7 

91 

•250869,31 

39482.47 

77463,56 

a 

9? 

-226325,44 

33599.91 

7950® . 38 

9 

93 

-200403.44 

27107.98 

63461.50 
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***************** VARI ANCF -CDVARI ANCE MATRICES ******************** 

THE vARlAMCf -CflV*Rl*NCE MATRICES FOR MODEL COORDINATES 


POINT 

7 

0 , 409600E 

02 

0,0 


0,0 




0.0 


0, 127690F 

03 

0,0 




0.0 


0.0 


0,28561 OF 

03 

POINT 

8 

0.504100F 

0? 

0.0 


0,0 




0,0 


0.386090F 

03 

0.0 




0,0 


0.0 


0,97344 OF 

03 

POINT 

9 

0.396900E 

0? 

0.0 


0.0 




0.0 


0. 106090F 

03 

0,0 




0,0 


0.0 


0.213160E 03 

POINT 

49 

0.409600F 

0? 

0.0 


0,0 




0,0 


0.280900F 

0? 

0,0 




0.0 


0.0 


0.268960E 

03 

POINT 

50 

0.640000E 

0? 

0.0 


0,0 




0,0 


0.435600T 

0? 

0.0 




0.0 


0.0 


0.B64360E 

03 

POINT 

51 

0.396900E 

u ? 

0.0 


0.0 




0.0 


0.291600F 

0? 

0.0 




0,0 


0.0 


0 , 1 82250E 

03 

POINT 

91 

0.739600E 

02 

0.0 


0,0 




0,0 


0.108160F 

03 

0,0 




0,0 


0.0 


0.412069F 

03 

POINT 

9? 

Oo 129980F 

03 

0,0 


0,0 




0,0 


0.243360F 

03 

0,0 




0.0 


0.0 


0. 1004B9F 

04 

POINT 

93 

0.577600E 

0 ? 

0.0 


0,0 




0,0 


0.&M900F 

0? 

0.0 




0,0 


0.0 


O.313290F 

03 


The variant -covariance matrices for ground coordinates 


0.0 0,0 

o.iooooor-i<? o.o 

0,0 0,100000F*19 


0,1 OOOOOF-19 

0,0 

0,0 
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WEIGHTS used in solution 


THE WEIGHT MATRIX FOR THE TRANSFORMATION PARAMETERS 


o,ioor-o5 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0.0 

0 , i 00 E *05 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0 . 100 E -05 

0,0 

0,0 

0,0 

0*0 

0,0 

0,0 

0,0 

0 . 100 E -05 

0,0 

0,0 

0,0 

0,0 

0.0 

0,0 

0,0 

0 . 1 OOE -05 

0,0 

0,0 

0,0 

0.0 

0,0 

0,0 

0,0 

0 . 100 E - L 5 

0,0 

0.0 

0,0 

0,0 

0,0 

0,0 

0.0 

0 . 100 F -05 


THF GROUND COORDINATE 

WflUE-T 

MATRIX 


0.1000F 21 

0,0 


0,0 

0.0 

0.100CF 

21 

0,0 

0,0 

0,0 


O.lOOOE ?1 


The WEIGHT MATRICES F°R The COORDINATES Of The mOODEL POINTS ( CONTROL AND UNk 


POINT 

7 

0.2441E-01 

0.0 

0.0 



0,0 

0.TP31E-02 

0,0 



0,0 

o.o 

0.3501E-02 

POINT 

6 

0.1984E-01 

o,o 

0.0 



0,0 

0.2577E-O? 

0,0 



0,0 

0,0 

0.1027E-02 

POINT 

9 

0.2520F-01 

0,0 

0.0 



0,0 

0.9A26E-0? 

0,0 



0,0 

O.C 

0 , 4691 £"02 

POINT 

49 

0 t 244lE-01 

o.c 

0,0 



0,0 

0.3560E-01 

0.0 



0,0 

0,0 

0.3718E-02 

POINT 

50 

0 # 1563E-01 

0,0 

0.0 


0.0 

0.2296E-01 

0,0 



0,0 

0.0 

0.115FE-02 



100 


POINT 

51 

0.2520F-01 

o.r 

0,0 



0,0 

0.0 

0.3A29E-01 

0.0 

0.0 

0.54P7E-02 

POINT 

91 

0,13520-01 

0,0 

0.0 

0.92«ftE“0? 

0,0 

o.o 



0,0 

0,0 

0.2427E-02 

POINT 

9? 

0.7695E-02 

0,0 

o.c 

0.4109E-0? 

0,0 

0,0 



0,0 

0.0 

0.9951E-03 

POINT 

93 

0,1731F>01 

0,0 

0,0 

0.0 

0.1 156E-01 
o.c 

0.0 

0,0 

0.3192E-02 
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I 1 E PAT I ON NUMBFR 1 



* SIGMA ZERO* 0 , 1 366E 00 * 


* * 


transformation parameters 


OLD VALUE CORRECTION NEN yALLE 


SCALE FACTOR 
TRANSLATION IN X 
TRANSLATION in y 

transi ation IN I 
angle omega 
angle phi 

ANGLE CAPA 


1 .0*620 
39.4*5 
47.348 
27,776 
0-11-1?, 32 
0-39-59.26 
0-36-36 ,79 


•0,4 6 ?1e- 01 
-0.3713F 02 
-0.45?9f 02 
-0.2A75F 02 
0-11-59,78 
0 0-30.43 

0 -4-15,99 


0 , 99999 
2.354 
2.059 
3.022 
0-23-12,10 
0“40“29.69 

0-40-52, 7fc 


ROTATIONS In RADIANS 


angle omega 
angle phi 

ANGLF CApA 


-0.00 J 2 ?95 0 

-0.01163193 

-0,01065033 


-0 . 34?0f" 02 
-0* 1475f -03 
-0. 12AU-02 


- 0*00674911 
-0«01 177945 
•0.01189142 


DEVIATION 

0 • 84r 3e -05 

O, 37?1 f 01 

P, ?569f oi 

0* 8 1 59f ol 

0 0 4, si 

P 0 7,87 

0 0 1,34 


0 . 2236^-04 

O. 38i 4f-o4 

P. 65o9f-o5 



* * * * * 
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1 1 f PM 1 ON NUMRFR 2 



* SIGMA 2PK0= 0 » 1 6 a 8E -0 1 * 



transformation PARAMETERS 


OLD VALUE CORRECTION NEN VALLE PF v I A T I ON 


SCALE F AC 1 OR 


0.9R999 

0 . 

67S2E-05 

1 ,00000 

O.10J7F-05 

TRANSLATION IN 

X 

2.35* 

-0. 

2576F 01 

-0,171 

o,as90E 00 

translation in 

Y 

2.0*9 

-0, 

1879f 01 

0.17V 

P.3242F 00 

TRANSLATION In 

2 

3 » 0?2 

-0. 

3 1 A 1 p 01 

■0,13V 

0.1 0?9f 01 

angle omega 


0-23-J ?. 1 0 

0 

0*33, 1 A 

0-23-45,2* 

0 0 0.58 

anglf phi 


0-*0"?9,69 

0 

0 *2.77 

0“*0"3?.Afc 

0 0 0.99 

angle CApA 


0-*0-82,78 

0 

0-11,6* 

0-* 1 "*,*2 

0 0 0.17 

POTATIONS IN RADIANS 





angle OMEGA 


-0, 0067*91 1 

-0. 1606F-03 

-0.00690976 

0.2870E-05 

ANGLF PHI 


-0,01 1 77945 

-0, 

1 3A *F-0* 

-0.0117928V 

p.abioe-o* 

ANGLF CAPA 


-0.011P9142 

•Of 

56A3F-0* 

-0.0119*78* 

P.81 80F-0* 
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FINAL OUTPUT 


The SOLUTION CONVFRGFS AFTER iteration number 3 


******************************** 
* * 

* SIGMA ZERO* 0.1638E-01 * 

* * 
******************************** 


TRANSFORMATION PARAMETERS 


OLD VALUE 


CORRECTION FINAL VALUF DEVIATION 


scale factor 

TRANSLATION IN x 
TRANSLATION IN Y 

translation in j 
angle omega 

ANGLE PHI 
ANGLE CApA 


1 tCCOOO 
-0.171 
0,*79 
“0.139 
0-23*45,24 
0 -Ao"3?.A6 
0 " A | - A . A 2 


0.80MF-06 

-0.78A0F-Q1 
-0 t ?03 6 p 00 
-0.38P6F-01 
0 0 *0.05 

0 0 0.11 
0 0 - o.ifl 


1.00000 
-0.25L 
- 0 . 0 ?* 
-0,177 
0*23*45.29 
0*40*32.34 
0” A 1 “4.6c 


O.lOllE-05 
O.A663E 00 
0,3??3e 00 
0.1 0?3f 01 
P 0 0.58 

P 0 0.99 

P 0 0.17 


the vApi ance-covarUncf mApTIx pf the transformation 


0 » 10?3 E -1 1-0.2A8i e -d 6-0» l979f-o8-o. ?6oOF*08-U.4o96p-l 3 0 » 4 372fl2 p. 14 38^12 
-0»?A81r*08 0 . 2 1 74 p 00 0*7l56f*0l 0»3997p 00 0. 7 o97 p- 06 “ 0. 1 98o£* 05 o»1033e"07 

-0*1 9 80f- 08 0. 7 1 56p*0l 0«1 o3 9 f 0(1 0.?118f 00 0 . 7 7 30p- 06-0, 9265 p* 06 0.17lBp-06 

-0«26o6p-o8 0 . 3997 f 00 0»?11 8 F 00 0»1046p 01 0*1961 E" 05-0.4634 t* 05 0»l823F-06 

- p .4o96f-13 0 . 7 097 f -o 6 0* 7730p*0« 0.1981 e* 05 0. 7«58 p- 1 1- 0. 8 556 p* 1 1 0»4640E“12 

0.4 372f-1?-0. 1 98 of* 05- 0 • 9265 f 06*0. 48 34 f- 05-0. 8558 E-n 0. 2L&7 t* 1 0* 0. 7759 p- 1 2 
0« 1 4 38 p- 1 2 0 . 1033f* 07 0*l7l8 f -o6 0. 1823 f-06 0. 4 64 OE* 1 2- 0 . 77 59 p- J 2 0.66llp*l2 



LIST OF ADJLSTF.D GROU^O COORDINATES 


NllMPFR X 

Y 

Z 

7 

“207558.638 

-81965.938 

P2P01 .31 3 

8 

- 165975, boo 

-86656, 3^5 

86663.375 

9 

-183591.375 

-91797.563 

86196.063 

69 

“228680, 000 

-2426o,516 

62161.063 

50 

“205500.813 

-29619.055 

86266,688 

51 

-181393.375 

-35455.576 

89601.500 

91 

“250889.313 

39682,669 

7 7663,563 

9? 

-226325.A3 8 

33599 . 9 H 

79504,375 

93 

-200903, A38 

27167.986 

83461,500 


SIG(X) 

S I (i f Y ) 

STG(Z) 

0. 1638f;-11 

0. 1638? “i 1 

0. 1638E-1 1 

0. 1 6 38 £• 1 1 

0» 1 6 30 F - 1 1 

0. 16 38 E* 1 1 

0. 1636F-11 

0. 1638 F" 1 1 

0. 16 36 E- 1 1 

0. I638f-n 

0 . i 6 SB F" ! 1 

0. 1638E-1 1 

0» 1 638p“ i l 

0 » 1 6 38 F- ) 1 

0.1838E-11 

0. I638p-ll 

0. 1638 E- 1 1 

0. 1 6 38 E* 1 1 

0. 1 638^“ 1 1 

0 * 1 6 3 8 F * 1 1 

0. 1*38 F-1 1 

0.1638E-11 

0, 1 63 8F- 1 1 

0.1638E-11 

0.1636^-11 

(m1*38E"11 

0. 163BE-11 
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